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Abstract 

We study two widely used algorithms for the Potts model on rectangular subsets of 
the hypercubic lattice Z d - heat bath dynamics and the Swendsen-Wang algorithm - 
and prove that, under certain circumstances, the mixing in these algorithms is torpid 
or slow. In particular, we show that for heat bath dynamics throughout the region of 
phase coexistence, and for the Swendsen-Wang algorithm at the transition point, the 
mixing time in a box of side length L with periodic boundary conditions has upper 
and lower bounds which are exponential in L^ 1 . This work provides the first upper 
bound of this form for the Swendsen-Wang algorithm, and gives lower bounds for both 
algorithms which significantly improve the previous lower bounds that were exponential 
in L/(logL) 2 . 

1 Introduction 

Convergence to equilibrium of heat bath dynamics and other dynamics for several lattice 
spin models of statistical mechanics has been of significant interest, for well over a decade, 
in probability theory, statistical physics, combinatorics and theoretical computer science. 
While the excellent monograph [31] provides a testament to this, many exciting new results 
and new techniques have since been developed. Fine examples of this development include, 
on the fast mixing front, results for Glauber dynamics on trees [5], for Swendsen-Wang 
algorithm on various classes of graphs [T7] , for a simple random walk on the super critical 
percolation cluster (U [23] . On the slow mixing side, results on the Swendsen-Wang for the 
Potts model on the complete graph [21], the heat bath algorithm for the Ising model at low 
temperature [30], and on quasi-local algorithms for the hardcore lattice gas model at low 
temperature [8] form similarly interesting and technically challenging examples. 

In this paper, we study two Monte Carlo Markov chains (MCMC), heat bath dynamics 
and the empirically more rapid Swendsen-Wang algorithm, for the g-state Potts model. Our 
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work is a continuation of the work begun some time ago in collaboration with several other 
authors [8], where we obtained weaker bounds than those we establish here. 

The point of our previous work was to relate the mixing times of MCMC in several models, 
including the Potts model, to the phase structure of the underlying equilibrium models. The 
Potts model is known to undergo a phase transition from a so-called disordered phase with 
a unique equilibrium state to an ordered phase with the coexistence of multiple equilibrium 
states. In our previous work we showed that, for the g-state Potts model on rectangular 
subsets of the hypercubic lattice Z d with periodic boundary conditions, heat bath dynamics 
is slow or torpid throughout the region of phase coexistence, while the Swendsen-Wang 
algorithm is torpid at the transition point, provided that q is large enough. There the lower 
bounds on the mixing time in a box of side length L with periodic boundary conditions were 
exponential in L/(logL) 2 . In this paper, we show that the mixing is even slower, obtaining 
essentially optimal results: both lower and upper bounds on the mixing time which are 
exponential in L . 

Slowness of the Swenden-Wang algorithm for the Potts model at the transition point was 
proved first on the complete graph [23]. This result initially came as a surprise to many 
physicists who had tacitly assumed that the algorithm was fast at all temperatures. Our 
previous work [8] was the first to establish such a result on subsets of the hypercubic lattice, a 
case which is both more physically relevant and technically much more challenging than the 
complete graph. To overcome these difficulties, we used some deep results from mathematical 
physics, which we now extend. In particular, our work brings to bear and extends, statistical 
physics expansion techniques for the problem of controlling the number of cutsets in graphical 
expansions of these models. Specifically, we use the so-called Pirogov-Sinai theory [36] from 
the statistical physics literature, in the form adapted to the Potts model by Borgs, Kotecky 
and Miracle-Sole ([ID], [II])- We also use the isoperimetric inequalities of Bollobas and 
Leader [6], as well as a large deviations technique borrowed from [7]. 

For Markov chains that change the value of only a bounded number of spins - such as 
the heat bath algorithm, it is easy to obtain upper bounds exponential in L d ~ x using either 
refined canonical path arguments as in [5] or recursive bounds on the Dirichlet form as in |14j . 
see [33] for a review. However, for the Swendsen-Wang algorithm, which is highly non-local 
in the spin representation, such an upper bound is not obvious. However, it turns out that 
a refinement of the bounds on the Dirichlet form can be used to obtain the desired upper 
bound. More generally, introducing a new graph parameter which we call the "decomposition 
width," we derive upper bounds for Swendsen-Wang on arbitrary graphs, which as a special 
case proves that Swendsen-Wang on trees is polynomial in the number of vertices for all 
temperatures. 

The real challenge is to obtain a lower bound which is exponential in L^" 1 , significantly 
improving the e cL//log L lower bound of [8]. While the previous bound required that we 
consider only contours which can be embedded into Z d , this optimal lower bound requires 
that we deal explicitly with the topology of the torus. In particular, we must distinguish 
between surfaces with vanishing and non-vanishing winding numbers, which we call contours 
and interfaces respectively. Moreover, in order to deal optimally with the contours, we need 
to define an appropriate notion of exterior and interior which allows us to establish a partial 
order on the set of contours. This in turn is used to develop the appropriate Pirogov-Sinai 
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theory on the torus. 

In order to state our results precisely, we need a few definitions. Let G = (V, E) be a 
finite graph and let (3 > 0. For a positive integer q, let [q] = {1,2,..., q}. The Gibbs measure 
of the (ferromagnetic) g-state Potts model on G at inverse temperature (3 is a measure on 
[q] v with density 

HG{°) = —y , (1-1) 

where cr £ [q] v is a spin configuration. Here 

Hq{<t)= (1-S(a x ,a y )) (1.2) 

xyeE 

is the Hamiltonian, and the normalization factor 

Z G = ^ HG(a) (1-3) 

is the partition function. (In the above, 5 denotes the Kronecker delta function.) 

For a finite A C Z d we define the measure /xa,i with the "1-boundary condition" by 
setting all the spins at the external boundary of A to 1. Explicitly, for cr £ {1, 2, ... , q} A , let 

H Atl (tr) = H GlA] (tr) + ^ (1-5(^,1)), 

y£A c 
\a>— Jl|=l 

where G[A] is the induced subgraph of A and | • | is the l\ distance. 
Then is the probability measure with density 

Va,i(<t) = — , (1.4) 

where Za,i is the appropriate normalization factor. 

The infinite volume magnetization is defined as M(/3) = lim^oo M\ L (/3) where A^ = 
{l,...,L} d and 




By standard correlation inequalities, the limit lim^oo M\ L (f3) is known to exist and to be 
monotone nondecreasing in (3. The transition point is defined by 

Po = (3 (Z d ) = inf{/3 : M(f3) > 0}. 

Let d > 2. Then < (3q < oo, and the model has a unique (infinite-volume) Gibbs state 
for (3 < (3o, and at least q extremal translation-invariant Gibbs states for /3 > /3 PO- For q 
small enough (depending on d), the model is believed to have a unique Gibbs state at (3 , 
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whereas for q large enough (again depending on d), it is known [2?] 129] to have q + 1 extremal 
translation-invariant Gibbs states at (3 , with 

P = ±\ogq + O(q^ d ). (1.5) 

Next we define the mixing time of a finite Markov chain with state space Q. Let P 
denote the transition probability matrix of an irreducible Markov chain with the (unique) 
stationary measure /i. The (variational) mixing time of such a chain is defined as 

r = min jt : d(t) < ^-J , (1.6) 

where 

d(i) = maxmax pi(A) — / P t ((T,(r') 

In this paper, we will consider several Markov chains for the Potts model on the torus 

T L4 = (Z/LZ) d . 

The chains we consider are the heat bath (or Glauber) dynamics, and the presumably much 
faster SW algorithm; see Section [27T1 of the definition of these chains. We denote the mixing 
times of these algorithms for the g-state Potts model on the torus T^d by r HB = r HB (Tirf) 
and r sw = r sw (T ijd ). 

Theorem 1.1 There are universal constants ki,k 2 < oo such that, for (3 > 0, d > 2 and 
any positive integer L, the following bounds hold 

r HB (T L4 ) < e ^ +k ^ Ld - 1 (1.7) 
r sw (T L4 ) < e ^+W)L d -\ (1.8) 

In order to prove this, in Section [31 we introduce the "partition width" of a graph, a 
notion that may be of independent interest. As a corollary of our proof, we also obtain that 
the mixing time of SW on a tree with n vertices, maximum degree <i max , and depth O(logn) 
is bounded by n l+ °^ duiax -\ see Corollary 13.21 for the precise statement. This generalizes 
the result of [5], which gives polynomial mixing for the HB algorithm on trees, to the SW 
algorithm. 

Theorem 1.2 Let d > 2. Then there exists a constant k% = ks(d) > such that, for q and 
L sufficiently large, the following bounds hold: 

r HB (T L4 ) > e^" 1 for all (3 > (3 (Z d ) (1.9) 



- SW (T ijd ) > e k ^ L ^ for (3 = (3 (Z d ). (1.10) 

Very roughly speaking, the reason for the heat bath lower bound is that this algorithm 
cannot move quickly among (the finite analogs of) the q translation- invariant extremal states 
present for > /3q. On the other hand, the Swendsen-Wang algorithm can move quickly 
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among these q states, but cannot move quickly between these q states and the one additional 
translation- invariant extremal state present at /3 = /?o- Whenever the algorithm cannot move 
quickly between states, in order for the system to mix, it must pass through a configuration 
with a "separating surface" of size at least L d ~~ x between the relevant states. 

The organization of this paper is as follows. In Section [2} we define the algorithms and 
the necessary notions from the theory of MCMC. In Section [31 we introduce the notion of 
partition width and establish the upper bound on the mixing time. In order to obtain the 
corresponding lower bound, we need some preparation: in Section [4] we construct the contour 
representation of the model, while Section [5] establishes the required geometric properties of 
contours and interfaces. This section can perhaps be skipped on first reading. In the next 
section, we state the necessary bounds from Pirogov-Sinai (to be proved in the appendix), 
and use these to establish two key estimates needed for the main proof: a bound on the 
probability of interfaces, and a suitable large deviation bound. Using these bounds, we then 
prove our main result in Section 

The reader only interested in the main flavor of our proofs should perhaps start with Sec- 
tion [2131 where we explain the main proof strategy, and then immediately jump to Section 
glancing back at Section H] and Section 16.11 as necessary. 



There are several MCMC algorithms that are used to generate a random sample from the 
distribution corresponding to the ferromagnetic Potts model. The heat bath is perhaps the 
simplest such Markov chain. Its transitions are as follows: Choose a vertex at random, and 
modify the spin of that vertex by choosing from the distribution conditional on the spins of 
the other vertices remaining the same. In contrast, the Swendsen-Wang algorithm can alter 
the spins on many vertices in each iteration. 

Throughout this section G = (V, E) is a fixed finite graph. For a subgraph G = (V, E) of 
G, we denote the set of (connected) components of G by C(G) = C(V, E), and its cardinality 
by c(G) = c(V,E). Finally, for a spin configuration er G [q] v , let E(cr) be the set of 
"monochromatic edges" xy G E with a x = a y . 

Heat Bath: From a spin configuration er G [q] v , we construct a new configuration er' as 
follows: 

HB1 Choose v uniformly at random from V. 

HB2 Take a' w = a w , for all w G V \ {v}, and change a v to a' v with probability 



2 MCMC Preliminaries 



2.1 Algorithms for the Potts Model 




(<t' v \(Tv\{v}) 



idE v 
vw£E 



g 




k=l 



■wev 



5 



For future reference we denote the transition matrix of this chain by 



p ™( cr > cr/ ) = iyrE(^ G Kknw) II ^^'^ 

In practice, an alternative method, the Swendsen-Wang algorithm [32], is often preferred. 
Swendsen-Wang Algorithm: For cr G [q] . 

SW1 Let E(cr) C E be the set of monochromatic edges. Delete each edge of E(cr) inde- 
pendently with probability 1 — p, where p — 1 — e _/3 . This gives a random subset 
A C E(cr). 

SW2 The graph (V, A) consists of connected components. For each component, choose a 
color (spin) k uniformly at random from [q], and for all vertices v within that compo- 
nent, set a' v = k. 

Again for future reference, we denote the transition matrix of this chain by 

AcE(cr) CeC(V,A) ^ k=l veV(C) 

The Swendsen-Wang algorithm was motivated by the equivalence of the ferromagnetic 
g-state Potts model and the random cluster model of Fortuin and Kasteleyn [22], which we 
now describe. Fortuin and Kasteleyn realized that the Potts model partition function Zq 
and expectations with respect to the measure \iq can be rewritten in terms of a weighted 
graph model on spanning subgraphs (V, A) C G with weights 

u G (A) = — p \ A \(l-p)\ E \ A \ q < v > A \ (2.1) 
Zg 

The relationship between the two models is elucidated in a paper by Edwards and Sokal 
[2T] . The Potts and random cluster models are defined on a joint probability space [q] v x 2 E . 
The joint probability vr(cr, A) is defined by 

*q(*,A) = -^ |A| (1 -p) 1 ^ 1 J] 5fa,*y). (2.2) 

By summing over a or A we see that the marginal distributions are vq or \iq respectively. 

A step er — > cr' of the Swendsen-Wang algorithm can be seen as (i) choose a random A' 
according to ^ci^i ■)// i G'( cr ) an d then (ii) choose a random cr' according to 7r(-, A')/i/g(A'). 
After Step SW1, we say that we are in the random cluster representation of the chain. 
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2.2 Mixing Time and Related Quantities 

Throughout this section, let P be the transition matrix of an irreducible Markov chain on a fi- 
nite state space Q, let \i be the stationary distribution of P, i.e., /i(cr') = J2a-en A t ( cr )-P( cr 5 °"')> 
and let /i m ; n = muio-g^ At(cr). We denote the mixing time, defined in ( 11 .61) . by r(P), and define 
the inverse gap (or eigentime) f(P) as 

W =T^ (2 ' 3) 

where the supremum is over all real valued functions g on Q with Var g > 0. Here as usual, 
Var g = Var^ # = - ^ (#(er) - #(er')) 2 ji{a) fi(cr') , 

and 

£(0,0) = £p(g,g) = 2^(0(°O -^(°" / ))V( cr )^( cr > cr/ )- 

<T,(T' 

If P is reversible, f (P) is just (1 — /^(P)) -1 , where /^(P) is the second largest eigenvalue of 
P. 

It is well known that the inverse gap can be bounded above in terms of the mixing 
times; if the chain is lazy, i.e., if the minimal self-loop probability muio- P(er, u) is uniformly 
bounded from below, a bound in the opposite direction is also not very hard to prove, see, 
e.g., [2]. However, the SW chain is not lazy. Instead of the standard upper bound on r(P) 
in terms of r(P), we therefore use a bound from [35]. For reversible chains, this bound gives 



r(P) <r(P 2 ) log(— ), (2.4) 



where, as usual, P 2 (cr, a') = Xl<r"ec P ( cr ' ( T ")P( (T ", &'), denotes the transition matrix of the 
two-step chain. 

We will also need an identity for the mixing time of a product chain. Let fli, Sl 2 be finite 
sets, and let P\,P 2 be the transition matrices of two irreducible Markov chains on f2 x and 
O2 with stationary distributions \x\ and /12 respectively. Let Q± x Q 2 be the set of all pairs 
a = (cr^, <x( 2 )) with E f2i and cr^> E Then the product chain is defined as the 
chain with the transition matrix 

{P^ P 2 ){<r^) = P l {a^\^)P 2 {^\^). (2.5) 

Let 

(/ii X fi 2 ){cr) = /i 1 (cr( 1 ))/i 2 (cr( 2 )). 

If both Pi and P2 have non- negative eigenvalues, then Pi x P 2 has non- negative eigenvalues 
and (3 2 (Pi x P2) = rnax{/32(Pi), /92(P2)}- Using this fact, one immediately shows that Pi x P 2 
is irreducible with stationary distribution fii x /i 2 and obeys the bound 

f(P x P 2 ) = max{f(Pi),f(P 2 )}. (2.6) 
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For our lower bounds, we use the notion of conductance and its relation to the mixing 
time. Setting 

Q(S,S')=J2 E^)^'^)' 

a-eScr'eS' 

the conductance of a set S C f2 can be defined as 

Finally the conductance of a Markov chain with the transition matrix P is 

$(P) = min $ 5 . (2.8) 

S:0<A»(S)<1 

The mixing time can be easily bounded from below in terms of the conductance, see 

s ^W)- (2 - 9) 

2.3 Proof Strategy 



In order to prove Theorem II. 1[ we will want to give an upper bound on the inverse gap f 
defined in (12.31) . To this end, it will be convenient to consider the SW algorithm on a general 
graph G. We then iteratively partition the set V of vertices of G into two sets V\ and V2, 
and show that the inverse gap of (Pj w ) 2 is bounded by the inverse gap of the product chain 
for SW on the two induced graphs G[Vi] and GfV^], times a factor which is exponential in 
the number of edges with one endpoint in V\ and one endpoint in V^. With the help of (12.61) 
and (12. 4ft . this allows us to bound the mixing time of SW by a number which is exponential 
in a quantity we call the partition width of the graph G. Applied to the torus Tl^, this gives 
a bound which is exponential in L d_1 , and applied to a tree, this will give a bound which 
is polynomial in the number of vertices. The bound for the HB algorithm follows a similar 
strategy. 

To prove Theorem 11.21 we w iU use the lower bound (12.91) . To this end, we will want to 
find a set of spin configurations S such that $5 is exponentially small in L d ~ x . Recalling 
that the Potts model at the transition point exhibits the coexistence of q ordered phases 
and one translation invariant phase, we will want to exploit the fact that the SW algorithm 
cannot transition easily between (the finite volume analogue of) the ordered phases and 
the disordered phase. To make this precise, we define S to be the set S = {cr : |_E(cr) | > 
(1 — a)dL d }, where a > is a small constant, say a = 1/3. Thus S consists of the 
configurations whose set of monochromatic edges form almost all of E. 

For large q the inverse transition temperature Po is large as well, implying that at (3 = (3 , 
the probability of deleting an edge in the first step of the SW algorithms is small; starting 
from a configuration cr G S it is therefore unlikely that after one step of the Markov chain, 
the new configuration cr' is such that the number of edges in E(er') is much smaller than 
(1 — a)dL d . (The probability that it is smaller than, say adL d , is actually exponentially small 
in L d ). But it is also unlikely that a configuration cr' £ S c has a number of monochromatic 
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edges which is larger than adL , since both requirements together imply that the number of 
edges lies between adL d and (1 — a)dL d . This corresponds neither to an ordered phase (which 
would have more than (1 — a)dL d monochromatic edges), nor to a disordered phase (which 
would have less than adL d monochromatic edges), and thus to a configuration which has low 
probability. Thus with high probability a configuration cr G S leads to a new configuration 
cr' which is in S as well, showing that <fis is small. 

To make this quantitative, we will have to show that the weight of the configurations in 
So = {cr' : adL d < \E(cr')\ < (1 — a)L d } is exponentially small in L d ~ l . To this end, we will 
first switch to the FK representation (12. ip . and then describe an edge configuration A C E 
in terms of a hypersurface separating regions with edges in A from regions with edges in 
E \ A. We then decompose this hypersurface into connected components, some of which to 
be called contours, and others to be called interfaces. While contours can have small or large 
size, interfaces will always have size at least L d ~ l . 

To prove our desired bounds on //(So), we will show that the probability of a configuration 
with an interface is exponentially small in the size of the interface, leaving us with the 
analysis of configurations without interfaces. These will come in two classes: configurations 
describing perturbations of A = E (we denote the set of these configurations by O rd), and 
configurations describing perturbations of A = (to be denoted by fidis)- Our last step then 
consists of a large deviations bound showing that with high probability, for configurations 
A G fiord, the common exterior of a set of contours has size at least (1 — |a)L d , and 
similarly for configurations A G fidis- This will imply that with high probability, the set of 
monochromatic edges of an ordered configuration cr has at least (1 — a)dL d edges, while that 
of a disordered configuration has at most adL d edges. Put together, these estimates give the 
desired bound for //(So). 

Our approach differs in several aspects from the approach taken in [8], which led to a 
conductance bound that was exponentially small in L/(log 2 L). 

First, the bounds in [8] relied on a combination of Pirogov- Sinai theory with the finite-size 
scaling theory developed in [10] and [11], where contours were by definition objects that could 
be embedded in M. d . This allowed for an immediate application of standard Pirogov-Sinai 
results, but produced error terms that were only exponential in L, which is not strong enough 
for our current purpose. To avoid this problem, we define contours in a purely topological 
manner, by requiring that their Z 2 winding number with respect to the torus is equal to 
zero. This has the advantage that the objects which cannot be classified as contours (we call 
them interfaces) must have size at least L^" 1 , since all surfaces of smaller size have winding 
number zero. The price we have to pay is that the set of contours now contains objects like 
long tubes winding around the torus which cannot be embedded into M. d , preventing us from 
applying the standard Pirogov-Sinai machinery as more or less a black box. 

Instead, we will show that Pirogov-Sinai theory does not really rely on the topology of 
R d , but rather on the implied structure of partial order on contours. More precisely, it relies 
on the fact that for any set of pairwise non-overlapping contours and interfaces, this partial 
order leads to a Hasse diagram that is a forest - this is expressed in Lemma 15.5} see also 
Lemma [4.31 and Definition 14.41 

Second, we will use a large deviation bound obtained by adding an artificial magnetic field 
to the contour model, see Section I6U1 for details. This turns out to be much more efficient 



9 



than the iterative, combinatorial arguments from [8], allowing us to improve a bound that 
is exponentially small in L/(logL) 2 to one which is exponentially small in L . 

3 Upper Bound on Mixing Time 

In order to prove our upper bound on the mixing time, it will be convenient to prove a more 
general theorem, involving a new notion, called the "partition width" of a graph, which we 
expect may be of independent interest. We need some notation. 

Given a graph G = (V,E), we define a hierarchical partition V of V by first dividing 
V into two non-empty subsets V\, V 2 , and then successively subdividing each set with more 
than one element into two further subsets until all sets contain only one element. Each 
such partition can be described by a rooted binary tree as follows: the vertices are subsets 
of V, with the root being V, the leaves being the singletons {x}, x e V. In addition, we 
have the constraint that for any vertex V{ with children Va, V^, we have Vi — Va U Vi 2 and 
Vn H Vi2 = 0- If Vj has children Va and V i2 , we define the weight w-piVi) of Vi as the number 
of edges between Va and V i2 in G; if V, is a leaf, we set its weight to zero. 

For x G V, we now define the separation cost sep v (x) of x as the sum of all vertex 
weights along the path from the root to x. The cost of a partition V is then defined as 
sep('P) = max^gy sep-p(x), and the partition width of G is defined as oj(G) = minp sep(P). 

Theorem 3.1 For any finite graph G = (V, E), the mixing time of the SW algorithm obeys 
the upper bound 

r(P| w ) < e 5 ^ G) (2+ \V\\og2 + /3\E\). 

Before proving the theorem, we state (and prove) the following corollary, which illustrates 
the usefulness of our notion of partition width. 

Corollary 3.2 Let A be a rectilinear subset of 7L d , let T^ ^ be the d-dimensional torus of 
side length L, and let T be a tree on n vertices with maximum degree d max and depth clogn. 
Then 

r(P| w ) < e 45/3|A(A)l (2 + (log 2 + d/3) | A|), (3.1) 
r(P| L W J < e 75 ^ 1 (2 + (log 2 + d(3)L d ), (3.2) 
r(P|W) < n 5c/w max ( 2 + Q Qg 2 + pj n j _ (3.3) 

Here A(A) is the volume of A divided by the minimal side length. 

Note that the second bound of corollary implies the bound fll.Sp in Theorem 11.11 The 
bound (11.71) of this theorem can either be proved by generalizing Theorem 13.11 to the heat 
bath algorithm (the proof is actually easier for this case), or by using the canonical path 
techniques of [5]. 

Proof of Corollary 13.21 Let G = (V, E) be an arbitrary finite graph, and let V = ViUV 2 be 
a decomposition of V into two disjoint subsets. Using the definition of the partition width, 
one easily verifies that 

w(G) < \E 12 \ +max{u(G[V 1 }),co(G[V 2 })}, (3.4) 



10 



where Eyi is the set of edges between V\ and V2. Using this bound, it is easy to verify by 
induction, that for a tree T of maximal degree d max and depth D, one has 

u(T) < d mauX D, 

which in turn gives the bound in the corollary for trees. 

We are thus left with proving upper bounds for the partition width for cubic subsets of 
Z d and the torus T^. We start with a rectilinear subset of side-lengths L\ > L% > ■ ■ ■ > Lj, 
which we denote by [L 1; . . . , LJ. To this end, we first note that 

w(G') < u(G) , 

whenever G' is a spanning subgraph of G, implying that oo([Li, . . . , L^]) is non-decreasing 
in the side-lengths Lj. Consider a set of sidelength L\ > • • • > L& with L\ > 2. Using the 
bound (13.41) . the monotonicity of oj([Li, . . . , Lj\) and the fact that [Lj/2] < 2Lj/3, whenever 
Lj > 2, we then bound 



([L!,...,L d ]) <L 2 ---L a 



< ^2 

< . . . 

< ( 



2\ 2 



+ 



3 



(3.5) 



< 3Li ■ • • + w 



r— 1 
1 2 1 



where k is the smallest z such that Li > 2. Using this bound, it is now easy to prove by 
induction that 

lo([Lx, . . .,L d ]) < 9A([L u ...,L d ]), 

implying the desired bound for the inverse gap on rectilinear sets. 

Next, we would like to bound the partition width of the torus using the just estab- 
lished bound for rectilinear sets. To this end we successively cut the torus in the d different 
coordinate directions, proceeding as in the proof above. Here, however, since we have a torus 
rather than a rectilinear set, we need two cuts rather than one cut in each direction to obtain 
two components. Keeping this in mind, we get 



Co' 



(T L , d ) < 2(1 



2\2 

3- 



3 



L x ---L 



d-l 



+ 00 

< QL d ~ l 



r 2 i> r 2 i>---> r 2 i 



9L 



d-l 



15L 



d-i 



which implies the desired bound on the mixing time. 

The proof of Theorem 13.11 is based on the following lemma. 



(3.6) 



□ 



11 



Lemma 3.3 Let G = (V,E), let P G = P| w and G = (V,E ), where E is an arbitrary 
subset of E. Then 

f{P 2 G)<f{P 2 G ) 

Proof. Recall that a single transition of the SW dynamics consists of two steps. Given a 
current Potts configuration cr, the first step identifies connected components of color classes 
and performs random edge deletion with probability e _/3 independently for each monochro- 
matic edge. We denote the resulting configuration by E' . The second step assigns colors 
independently for each new connected component (cluster) in E', resulting in a new Potts 
configuration cr'. 

Let E(cr) = {xy G E : a x = o~ y }. Let G and Go be as in the statement of the lemma, 
and let E l =E\E . Then 



Pair,*') = £ (i-e-o is v^)^ n ~n 

E'CEO) CeC(V,E>) " x,yeC 

> e -m\ j2 {i- e ^f\^ E ^ E ^ E '\ n iu s ^< 

E'CE(tr)\Ei 



Q 

E'CE(<r)\Ei CdC(V,E') * x,y£C 



implying that 

P G {a^')>e-^P Go {rr^'). (3.7) 
Next we observe that, by the definition of H G , we have 

e -/9fr Go (°-) e -/3|-El| < e -/3H G (cr) < e -(5H Go (cr) _ 

implying that 

Combined with ( 13. 7p and the definition of variance and the Dirichlet form, this proves that 

Var MG 9 < g^Var^ g < Var^ g 



and hence 



£p*(g,g) ~ e-mEA £ - e ^ g y 



G 

□ 

Having established Lemma 13.31 we are now ready to prove Theorem 13.11 
Proof of Theorem 13.11 Given a graph G = (V, E) and a decomposition of V into two 
disjoint subsets V\ and V 2 , let E i2 be the set of edges in E that join V\ and V 2 , and let 
E = E\ E 12 . Let d = G\Vi], G 2 = G[V 2 ] and G = (V, E ). Observing that no edge in G 
joins V\ and V 2 , we clearly have P Go = P Gl x P G2 and thus P Go = P Gi x Pq Combining 
Lemma [3.31 with the identity (12.61) for P G and Pq , we thus have 

f{P G ) <max{f(P^),f(P| 2 )} e 5 »l. 
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Applying this bound recursively for the decompositions in a hierarchical partition V of G, 
we obtain that 

f{Ph)<e^Y[f{Ph [{x}] ). 

x£V 

Since the inverse gap for the SW algorithms on a single site is one, we get t(Pq) < e 5 ^ G \ 
Combined with the bound (12.4j) . this proves the theorem. □ 



4 Contour Representation 

In this section, we derive a representation for the Potts model in terms of contours and 
interfaces. This representation is a modified version of the representation of [11]. We use 
T = to denote the d-dimensional discrete torus of sidelength L, with vertex set V = 
Vl,cL = (Z/LZ) d and edge set E = El^, and fl to denote the SW configuration space 2 E , i.e., 
the set of all edge configurations Ad E. 

We start from the random cluster representation (12. ip . Given A C E, let V(A) = 
^{x,y}eA{x, y}, and let G(A) = (V(A),A). Recalling the definition of c(A) as the number 
of connected components of the graph G(A) = (V,A), we introduce c(A) as the number of 
connected components of G(A). Observe that c(A) = c(A) + \ V \ V(A)\ and 

2\A\ = 2d\V(A)\ - -2\8 2 A\, 

where 

S k A = {{x, y} E E\A; \{x, y} fl V(A) \ = k} k = 1, 2. 

Using the notation SA = {{x,y} e E \ A; |{x,y}n V(A)\ > 0} and \\6A\\ = \SiA\ +2\5 2 A\, 
we rewrite the weight w(A) = ZtVt{A) of a configuration A in (12. ip as 

= q c(A) e -e di3 \V\V(A)\ e -e old \V(A)\ e -K\\SA\\ 

where 

e dis = -log((l -p) d q) =d(3- logg, (4.2) 
e ord = -rflogp = -dlog(l - e~P) = 0(e-P), (4.3) 

, = -\ log(^) = \ log(e^ " l) = f + 0{e-l>). (4.4) 

Throughout this section we will assume that (3 > log 2, so that, in particular, k > 0. 

The representation (14. ip already shows that there are three regions of interest: for /3 
^logg, the dominant configurations are those with most vertices belonging to V \ V(A), 
i.e., most vertices are isolated, corresponding to a "disordered high-temperature phase;" 
for (3 3> \ 1°S 9' ^ ne dominant configurations have most vertices in the "ordered region" 
V(A) corresponding to an "ordered low-temperature phase;" and for j3 ~ ^logg and q 
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(4.1) 



(and hence k) large, the dominant configurations fall into two classes, one with mostly 
isolated vertices and the other with most vertices in V(A). As we will see, even the SW 
algorithm has difficulties transitioning between these two classes, which leads to slow mixing 
at/3 = ilogg + O(g- 1 ^). 

We will decompose 5 A into several pieces called "interfaces" and "contours". More 
precisely, we will first "fatten" the set A into a subset V(A) of the continuum torus V = 
(R/ (LZ)) d and then decompose the boundary of V(A) into components, some of which will 
be called interfaces, and some of which will be called contours. 

For an edge e = {x, y} G E, let e be the set of points in V that lie on the line between 
x and y. Given A, we call a closed k- dimensional unit hypercube c C V with corners in 
V occupied if all edges e with e C c are in A. We then define the set V(A) C V as the 
1/4- neighborhood of the union of all occupied fc-dimensional hypercubes, k = 1, . . . , d, i.e., 
V(A) = {x G V : 3 c occupied, such that dist(x, c) < 1/4}, where dist(x, y) is the ioo- 
distance between the two points x and y in the torus V and dist(x, c) = infj, 6c dist(x, y). 
Note that the set V(A) is aunion of cubes of side-length 1/2 with centers in V\/2 = 'U£) d , 
and that the set V(A) defined at the beginning of this section is just the intersection of V(A) 
with the vertex set V of the discrete torus T. 

Alternatively, one can define the set V(A) by constructing its complement, the "dis- 
ordered region" V \ V(A) as follows: Let (E \ A)* be the set of d — 1 dimensional unit 
hypercubes dual to the edges in E\A. It is then easy to see that the union of the open 3/4- 
neighborhood of V \ V(A) and the open 1/4- neighborhood of (E \ A)* is just the disordered 
region V \ V(A) (see Lemma [5.11 in Section EJ). 

For i = 1, 2, . . . , d, let Lj be the fundamental loop = {y G V: y^ = 1 for all j ^ i}. If 
A is a union of cubes with diameter 1/2 and centers in V1/2 and 7 is a component of dA, 
then the winding vector N(7) G {0, l} d , with its zth component being equal to the number 
of intersections (mod 2) of 7 with Lj. 

Definition 4.1 Let A be a configuration in Q. The contours corresponding to the configu- 
ration A are defined as those connected components 7 of the boundary ofV(A) which have 
winding number zero, N(7) = 0; the remaining connected components of the boundary of 
V(A) are called the interfaces corresponding to A; the set of these interfaces is called the 
interface network corresponding to A. We denote the set of contours corresponding to A by 
T(A), and the interface network corresponding to A by S(A). 

Without reference to a configuration, we say that 7 is a contour if there exists a con- 
figuration A such that 7 G T(A), and similarly for an interface and an interface network. 
We define two contours (or two interfaces, or one interface and one contour) 7, 7' to be 
compatible, z/dist(7, 7') > 1/2. We define a contour 7 and an interface network S to be 
compatible if ^ is compatible with all interfaces in S. 

Note that the contours and interfaces corresponding to a configuration A G VL are ori- 
entable in the standard topological sense (see e.g. [12] or[3]); in fact, they are oriented, with 
an "ordered side" facing \(A), and a "disordered side", facing the complement of \(A). 
Thus our contours are "labeled contours" , with labels £ G {ord, dis} indicating which side is 
ordered, and which is disordered. 
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Note also that the contours and interfaces corresponding to a configuration A are pair- 
wise compatible. It is not true, however, that any set of pairwise compatible contours and 
interfaces correspond to a configuration A £ Q. In order to get a one-to-one correspondence, 
we define the notion of "matching labels." 

Definition 4.2 Let S be an interface network, and let T be a set of contours. We say that 
S U T is a set of matching contours and interfaces if the following conditions hold: 

(i) The contours and interfaces inTUS are pairwise compatible. 

(ii) The labels are matching in the sense that, for each component C of V \ U 7 eru.s'"}'' 
there exists a label £(C) £ {ord, dis} such that the ordered side of each contour (respectively, 
interface) faces an ordered component (i.e., a component with label £(C) = ord), and similarly 
for the disordered sides. 

For a set of matching contours and interfaces, we denote the union of the ordered com- 
ponents by V or d, and the union of the disordered components by Vdi S - 

With this definition, the set of contours and interfaces corresponding to a configuration 
A £ Q are clearly matching. It turns out (see Corollary 15.121 in Section [5]) that the converse 
is also true, namely that any set of matching contours and interfaces corresponds to exactly 
one configuration A £ Q. Using this fact, we rewrite the partition function Z as a sum over 
sets of matching contours and interfaces. 

To this end, we first note that the number of components c(A) is clearly equal to the 
number of components c(V(A)) of the continuum set \(A), and hence also to the number of 
components of the set V or( j. Note also that both V or( j and V^is are functions of the matching 
contours and interfaces (I\ <S). As a consequence, the weight (14. ip can be rewritten as a 
function of (T,S) according to: 

w(A) = q c ( v ( A )) e - e ^\ v \ v ( A )\ e - e oM\vnv(A)\ e -K\\dV(A)\\ 

= g c ( v ord) e -edis|v dis ny| e -e ord |v ord ny| "Q e -«||S|| e - K IMI ; (4.5) 

where ||<9V(v4)|| is the number of intersections oidY(A) with the continuum set E = U e gE e ' 
and similarly for \\S\\ and ||7||. Together with the already mentioned Corollary I5.12[ about 
the one-to-one correspondence between configurations and sets of matching contours and 
interfaces, this leads to the desired representation of the partition function Z = Zj- = 

1 = ^2 g c ( V °rf) e -edis|V dis nV| e -e ord |V ord ny| e -«||S[| g-relbi^ (4_g) 

5,r ses 7 er 

where the sum runs over sets of matching contours and interfaces. 

Next we define the interior and exterior of a contour. To this end, we need the following 
geometric lemma. Its proof is deferred to Section [5] 

Lemma 4.3 Let A be a configuration in Q, and let 7 be a contour of A. 

i) The set V \ 7 has exactly two components. 

ii) Let C and D be the two components of V \ 7, and let Si, S2 be two (not necessarily 
compatible) interfaces that are both compatible with 7. Then both Si and S2 lie either in C 
or D. 
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For the purpose of the next definition, it is convenient to define the size of a set W C V 
as the cardinality of W D V. 



Definition 4.4 Let 7 be a contour. If there exists an interface S (not necessarily corre- 
sponding to the same configuration) that is compatible with 7, we define the exterior, Ext 7, 
0/7 as the component o/V\7 that contains S; otherwise we choose the larger of the two com- 
ponents; finally, if both of these components have the same size, we choose that containing a 
distinguished point, x £ V. The interior is defined as the set Int7 = V \ (7 U Ext 7). 

Given a set of pairwise compatible contours T, we define a contour'-/ G T to be an external 
contour in T if there exists no contour 7' G T \ {7} such that Int7 C Int 7'. We also define 
the exterior of T as the set 



7er 

Finally, we say that 71,..., 7 n o,re mutually external if they are pairwise compatible and 
Int 7j D Int 7j = for all i ^ j . 

The next lemma states several properties of the set ExtT(v4) of a set of contours corre- 
sponding to a configuration A without interfaces. Its proof is again deferred to Section [51 

Lemma 4.5 Let A be a configuration with S(A) = 0, and let A ext be the set of edges with 
both endpoints in Extr(A). Then the following statements hold 

(i) Extr(A) is a connected subset of the continuum torus V. 

(ii) Ext T(A) is either a subset of the ordered region V(A) or the disordered region \\V(A). 
(Hi) //Extr(y4) C \{A), then (V(A ext ), A cxt ) is a connected subgraph ofT^^. 

5 The Geometry of Contours and Interfaces 
5.1 Elementary Topological Notions 

We start by reviewing some standard notions from algebraic topology. Let Ti/2 — {{^)/L%) d 
be the torus with two points connected by an edge if they differ by 1/2 (mod L) in one 
component. Its vertex and edge sets will be denoted by V1/2 and Eyi, respectively. 

We define /c-cells in V1/2 as the k- dimensional elementary cubes in V1/2, so that the points 
in V1/2 are 0-cells, the edges are 1-cells, etc. We also consider the dual V*/ 2 of V1/2, consisting 
of the barycenters of the <i-cells in V1/2, when considered as subsets of the continuum torus 
(R/LZ) d . As usual, given a fc-cell c in V1/2, we define its dual as the d — k cell c* in V*i 2 
that has the same barycenter as c, and similarly for the dual of a k-cel\ in Vf, 2 . Note that 



Given a /c-cell c, we define its boundary as the set of all (k — l)-cells that are subcubes 
of c (note that there are 2k such subcubes). More generally, for a set K of /c-cells, we define 
the (Z 2 -)boundary of K as the set of all (k — l)-cells which are in the boundaries of an odd 
number of cells in K. We denote this boundary by dK. The co-boundary of a k-ce\\ c is 
defined to be the set of (k + l)-cells which have c in their boundaries (there are 2(d — k) such 
(fc + l)-cells). 




(4.7) 



= c. 
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We often identify a fc-cell with the closed /c-dimensional continuum cube with corners 
being the vertices of the discrete cell. In this context, OK is identified with the corresponding 
continuum boundary. 

The (d — l)-cells in V*, 2 are called facets; the set of all such facets is denoted by F*, 2 . We 
say two facets are connected (or adjacent) if they share a (d — 2)-dimensional cell in their 
boundaries. 

A sequence of points xo,...,Xk G V1/2 is called a loop in T1/2 if £0 = Xk and {xi, xi + \} G 
E1/2 for all / G {0, k — 1}. Such a loop is called a fundamental loop in the i th lattice direction 
if k = 2L and all edges point in the i th lattice direction, and it is called an elementary loop 
if k — 4 and neither xq = X2 nor x\ = £3. Note that there are exactly 2d(2d — 2){2L) d 
elementary loops in T 1 / 2 . 

Consider now a set of edges X C Ei/ 2 and its dual X* = {e* : e G X}, and assume that 
X* is orientable. If dX* = 0, then we say that X* is an orientable closed surface, and we 
define the Z 2 -winding vector of X* as the vector N(X*) = (N h ...,N d ) G {0, l} d with JV* 
equal to the number of times X* intersects a fundamental loop in the i th lattice direction 
mod 2. 

5.2 Preliminaries 

Our first lemma summarizes several simple properties of the construction used in the defini- 
tion of contours. It involves both facets in V*, 2 , the objects dual to the 1-cells in Vyz, and 
(d— l)-dimensional unit hypercubes dual to the edges in E. While the first will be considered 
to be abstract objects in the sense of algebraic topology, the second will be considered to be 
closed hypercubes in V. We trust that this does not cause any confusion to the reader. 

We need some notation. Given a set D C E, let D* be the set of d — 1 dimensional unit 
hypercubes dual to the edges in D, let V-(D) be the set of all vertices x G V such that 
all edges {x, y} G E containing the vertex x lie in D (alternatively, this set can be defined 
as V-(D) = V \ V(E\ D)), and let Vdi S (-D) be the union of the open 3/4- neighborhood of 
V-(D) and the open 1/4- neighborhood of D*. Note that the set V-(D) and hence the set 
Vdis(-D) depends implicitly on E. Since E is fixed throughout, we suppress this dependence 
in our notation. 

Lemma 5.1 i) For ieO, the boundary d\(A) ofY(A) is regular in the sense that each 
(d — 2) -cell with corners in V*, 2 is shared by either zero or two facets in d\(A). 

ii) Let A G f2, and let C be a component of V \ dV(A). Then C is either a component 
ofV(A) or ofV\V(A). 

Hi) Let Ci, . . . , Cfc be the connected components of \(A), let Vi — V D Cj, and let Ai 
be the set of edges whose endpoints lie in V{. Then (Vi,Ai), . . . , (V^, A^) are the connected 
components of (V, A), and Cj = V(Aj). 

iv) Let D = E \ A and let Vd; s (-D) be as defined above. Then V \ V(A) = Vdi S (-D). 

v) Let C be a component of Vdi S (-D) ; and let Dq, be the set of edges in E whose midpoint 
lies in C. Then C = Vdi S (-Dc)- 

Proof, (i) Given a configuration Aefl, let V\/ 2 {A) be the intersection of V(A) with the 
vertex set Vi/ 2 of the discrete torus Ti/ 2 . The boundary of V(A) can then be rewritten as 
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the union of all facets that are dual to an edge e G E1/2 joining V\/2{A) to its complement 
in Vi/2- Using this fact, we easily prove the first statement of the lemma. 

Indeed, let e be a (d — 2)-cell with corners in V*, 2 , and let fi, f2, f3, ft be the four 
facets in the co-boundary of e. Then fi, f2, f3, are dual to edges £i£2, ^2^3, ^3^4, ^4^1, 
where the points £1, £2, £3, £4 comprise a closed path of length four in the torus T 1 / 2 . Let 
c(xi), . . . , c(x 4 ) be the four d-cells with centers xx, ... X4. Exactly one of these four cubes, 
say the cube c(xi), will be a cube whose center lies in the original vertex set V (recall that 
vertices in V have twice the spacing of those in 14/2)- Our convention of filling in only those 
hypercubes in the original torus T^d whose edges are all in A implies that either none or all 
or exactly one of the four cubes c(xi), . . . , c(xA lies in \(A). In the first two cases, none of 
the facets fi, f2, fs, Ja are in the boundary of \(A), and in the third case exactly two are in 
the boundary of V(A), which proves that each (d — 2)-cell with corners in V*, 2 is shared by 
either zero or two facets, as claimed. 

(ii) This is obvious. 

(iii) This now follows immediately from our fattening procedure for the ordered region, 
which respects the component structure of (V,A). 

(iv) We first prove that V dis (E \ A) C V \ V(A). Consider first an edge e G E \ A and 
its dual e*. We claim that all points with distance less than 1/4 from e* lie in V \ \(A). 
Since the set V(/l) increases if A increases, it is clearly enough to prove this statement for 
A = E \ {e}, in which case it follows immediately from the way we set up our fattening 
procedure for V(A). In a similar way, one proves that all points with distance less than 3/4 
from the vertices in V^(D) = V\ V(A) lie in V \ V(A). 

To prove that V\V(A) C V dis (E\A), let x G V\V(A), and let c C V be a d-dimensional 
unit cube with corners in V such that x G c. If x has distance less than 1/4 from the center 
of c, then at least one edge e G c(lV must lie in E \ A, since otherwise c would have been 
filled in our fattening procedure, contradicting x G V \ \(A). But the 1/4 neighborhood 
A/i/4(e*) of e* contains all points with distance less than 1/4 from the center of c, implying 
that x G Ai/4(e*) C Vdi S (£' \ A). If x has distance 1/4 or more from the center of c, it must 
have distance less than 1/4 from a d — 1 dimensional unit cube Ci in the boundary of c. If 
the projection, x±, of x onto Ci has distance less than 1/4 from the center of c 1; then one 
of the edges in Ci D V must lie in E \ A, which again implies x G A/1/4 (e*) C V dis (£' \ A). 
Continuing inductively, we are left with the case that x has distance at most 1/4 from one 
of the corners, y, of c. But this means that none of the edges containing y can lie in A, thus 
y G V \ V(A), and hence x G Af 3/4 (V \ V(A)) C V dis (£ \ A). 

(v) Noting that the midpoint of an edge lies in Vdi S (-E' \ A) if and only if its dual lies in 
Vdis{E \ A), this is an immediate consequence of our fattening procedure for the disordered 
region. □ 

Lemma 5.2 Let S\ and S2 be two interfaces with Si fl S2 = 0. Then N(Si) = N(S < 2). 

Proof. The interfaces S\ and S2 are closed orientable submanifolds of the torus (R/ U£) d . In 
the language of algebraic topology, the winding numbers N(S'i) and N(S , 2 ) are the Poincare 
duals of the submanifolds S\ and S2. The Poincare dual of the transverse intersection of 
two such submanifolds is then given by the wedge product of the Poincare duals of the 
submanifolds, see [12], Section 6. Since empty intersection is a special case of transverse 
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intersection, we conclude that the wedge product of N(Si) and N(5 < 2) must be zero. Let 
e*j be the unit vector whose j th coordinate is Sij. Recalling that e*j A e \ = and e*j A e} = 
— e} A e*j for all z, j, the condition N(Si) A N(S2) = is equivalent to the (2) conditions 
N i (S 1 )N j (S 2 ) - Nj(Si)Ni(S 2 ) = 0, which implies that N(Si) and N(S 2 ) are multiples of 
each other. Since both are different from 0, we conclude that N(Si) = N(S , 2). □ 

We close this section with an elementary lemma about "cutsets". It is best formulated 
in the context of a general connected graph G = (V,E). As usual, a subset E' C E is a 
cutset if (V, E \ E') is disconnected. It is called a minimal cutset if no proper subset of E' is 
a cutset. The following lemma is elementary; its proof is left to the reader. 

Lemma 5.3 Let G = (V, E) be a connected graph, and let W C V. If the edge-boundary of 
W is a minimal cutset then both the induced graph on W and V \ W are connected. 

5.3 Proofs of Lemma 14.31 and 14.51 

We start with the proof of Lemma 14.31 

Proof of Lemma 14.31 i) 7 has no boundary, is orientable, and has winding number zero. 
Therefore any closed path intersects 7 an even number of times, implying that 7 is the 
boundary of some open set C. Let D be the set D = V \ (C U 7). Then both C and D 
must be connected. Indeed, considering 7 as dual to a minimal cutset on the half-integer 
lattice, the connectedness of C and D follows immediately from the corresponding statement 
(Lemma 15. 31) for minimal cutsets. 

ii) The interfaces Si, S 2 are connected subsets of V. Since they do not intersect 7, each 
of them must lie in one of the connected components of V \ 7. We will have to prove that 
they both lie in the same component of V \ 7. 

Assume, by contradiction, that Si C C and S 2 C D. Since C fl D = 0, this implies 

51 n S2 = 0, which in turn, by Lemma [5.21 implies that N(Si) = N(S 2 ). Together with the 
fact that N(7) = we conclude that the winding number of 7 U Si U S 2 is zero, implying 
that any loop in V must intersect 7 U Si U S 2 an even number of times. 

Consider now a component Ni(Si) of the winding vector N(Si) that is not equal to 0, and 
let u be a fundamental loop in the i th direction, oriented in an arbitrary but fixed fashion. 
To reach our contradiction, we will modify u in such a way that it does not intersect 7 or 
S 2 , while intersecting Si an odd number of times. First, we note that the original loop u 
intersects Si an odd number of times. If it does not intersect 7, then it does not intersect 

5 2 either, since S 2 lies in D, while Si lies in C. If u; intersects 7, let x be one of the 
intersection points, and let y be the next intersection point (since N(7) = 0, the number 
of these intersection points must be even, so there must be such a y). Recalling that 7 is 
connected, we now replace the segment of u that joins x to y by a path in 7, and then 
deform this segment in such a way that it lies completely in C without intersecting Si (this 
is possible since dist(Si,7) > 1/2). Given that the original segment from x to y did not 
intersect Si since Si C C while the segment was a path in D, we have not changed the 
number of intersections with Si, so this number is still odd. Repeating this step until we 
have no intersections with 7, we remove all parts of u that lie outside of C, ending up with 
a path inside C that intersects Si an odd number of times, as desired. Note that the final u 
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may not be a lattice path, but nevertheless it gives the desired contradiction, since winding 
numbers are topological invariants in the continuum as well. □ 



Definition 5.4 Let 7 and 7' be two contours. We say that 7 and 7' are mutually external 
(and write 7 _L n f l ) if 7 is compatible with 7' and Int7 fl Int7' = 0, and we say that 7 lies 
inside 0/7' ('and wzte 7 < 7' J 2/7 zs compatible with 7' and Int7 C Int 7'. 

Lemma 5.5 Lei 7, 7' and 7" &e contours. 

(i) 1/7 and 7' are compatible, then exactly one of the following three holds: 7 < 7', 7' < 7, 



(it) Ifj± 7', iaen dist(Int 7, Int 7') > 1/2, and 2/7 < 7', then dist(Int7,Ext7') > 1/2. 

(mj // 7 < 7' and 7' < 7" inen 7 < 7". 
(iv)If 7 < 7' and 7' _L 7" iaen 7X7". 

Proof. (i) Let Int 7 and Ext 7 be the two components of V \ 7. Since 7' is connected 
and 7 fl 7' = 0, 7' must lie in one of the two components of V \ 7; therefore we have that 
either 7' C Int 7 or 7' C Ext 7. Since the same statement holds with the roles of 7 and 7' 
exchanged, we get that exactly one of the following four cases must hold: 



We claim that the last case is impossible. 

To prove this, we first show that (15.41) implies that Ext 7 fl Ext 7' 7^ 0. Recalling 
Definition 14.41 of Ext 7, note that Ext 7 is defined differently in several distinct cases. Let us 
first consider the case that there exists an interface S with dist(5, 7) > 1/2 and S C Ext 7. If 
dist(5, 7') < 1/2, then also dist(S, Int 7) < 1/2 by our assumption that 7' C Int 7. But this 
is not compatible with dist(S, 7) > 1/2 and S C Ext 7. Thus we have dist(S', 7') > 1/2. But 
if dist(5, 7') > 1/2, then S C Ext 7' as well, implying in particular that Ext 7 fl Ext 7' 7^ 0. 

Let us now consider the cases in the definitions of Ext 7 and Ext 7' such that there is 
no interface S compatible with either 7 or 7'. Consider the subcase that Ext 7 is defined 
by size, so that |Ext7 n V\ > |Int7 fl V\ and |Ext7' n V\ > n V\. This implies that 

|Ext7ny| > \V\/2 and |ExtynV| > \V\/2 which in turn implies that Ext 7 n Ext 7' ^ 0. 
The case where Ext 7' is defined by size is strictly analogous, so we are left with the subcase 
where both Ext 7 and Ext 7' are defined by containing the distinguished point Xq. Again, 
this implies that Ext 7 fl Ext 7' 7^ 0. 

Now we show that the condition Ext 7 fl Ext 7' 7^ rules out the case (15.41) . Let 
u G Ext 7 fl Ext 7', and let x G Int 7. Consider a path u from x to u, and let y be the 
last exist point from Int 7 along u. This implies y G 7, and by our assumption (15.41) . we 
therefore have y G Int 7'. But this implies there exists a point z G 7' C Int 7 after y, which 
is a contradiction. Thus case (15 .4p is impossible. 

Next we prove that the three remaining cases imply that Int 7' fl Int 7 = 0, Ext 7' fl 
Int 7 = and Int 7' fl Ext 7 = 0, respectively, showing that either 7 _L 7', 7 < 7' or 7' < 7, 
respectively. 



or 7 _L 7'. 



7 C Ext 7' and 7' C Ext 7 

7 C Int 7' and 7' C Ext 7 

7 C Ext 7' and 7' C Int 7 

7 C Int 7' and 7' C Int 7 



(5.1) 
(5.2) 
(5.3) 
(5.4) 
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We prove all three statements in one sweep, by setting A = Ext 7 and A' = Ext 7' in the 
first case, A = Ext 7 and A' = Int 7' in the second case, and A = Int 7 and A' = Ext 7' in 
the third case. Our assumption then reads 7 C A' and 7' C A, and our claim is B D -B' = 0, 
where I? = v4 c \ 7 and -B' = (A') c \ 7'. In a preliminary step, we prove that 7 C A' implies 
B\B' 7^ 0. Indeed, assume the contrary, i.e., B <Z B' . Taking the closure on both sides, 
this gives 7 U B C 7' U I?', and hence 7 C 7' U B' = (^4') C 5 a contradiction. Now we prove 
the main claim B D -B' = 0. Assume the contrary, that there exists an 2 £ I? D B' . From 
our preliminary claim, we also know that there exists a y £ B \ B' . Since B is connected, 
we conclude that there must be a path 00 G B from x to y. Let z be the first time this path 
exits B'. Then z £ 5 n dB' = B n 7' . Thus 5 n 7' 7^ 0, a contradiction. 

(ii) Both statements follow from the observation that if An A' = 0, then dist(cL4, cM') < 
dist(iM'). 

(iii) Assume that 7 < 7' and 7' < 7". Then Int 7 C Int 7' C Int 7", so we need 
only prove that 7 and 7" are compatible, i.e. that dist(7, 7') > 1/2. On the other hand, 
Int 7 C Int 7'. Taking the closure of both sides and using a trivial inclusion, this implies 
that 7 C 7 U Int 7 C 7' U Int 7'. Thus 

dist( 7 , 7") > dist(Y U Int 7', 7" U Ext 7") = dist(Int 7', Ext 7") > 1/2 

where we used (ii) in the last step. 

(iv) This is proved strictly analogously to the proof of (iii). □ 
The next lemma is an easy corollary of Lemma 15.51 

Lemma 5.6 Let T be a set of pairwise compatible contours. Then ExtT is a connected 
subset o/V. 

Proof. Let r ex t be the set of external contours in T. It follows immediately from the 
definition of ExtT and the last lemma that ExtT = Extr cxt . It is therefore enough to 
consider a set T of mutually external contours. We prove the statement by induction on the 
number of contours in T. The statement is trivial if T = 0. Assume the statement is proved 
for T = {7 1; . . . , 7n-i}. Adding an additional mutually external contour 7 n will not change 
the connectivity. Indeed, let x, y £ Ext (r U {7 n }), and let w be a path in Ext T that joins x 
to y. If u does not intersect j n U Int j n there is nothing to prove. Otherwise, let x' £ 7 n be 
the first entry point into 7„ U Int j n , and let y' £ 7„ be the last exit point from 7„ U Int j n . 
Since 7„ is connected, we can replace the path oj from x' to y' by a path in 7 n , leading to a 
path uj' joining x and y in 7„ U Ext (V U {7n}). By deforming u' this immediately leads to 
a path in Ext (r U {7 n }) 5 proving the lemma. □ 

We are now ready to prove Lemma 14.51 
Proof of Lemma 14. 51 (i) Since the contours corresponding to a configuration A are pairwise 
compatible, this statement follows immediately from the previous lemma. 

(ii) Using the fact that the boundary of V(A) is equal to the union over all contours in 
T(A), we conclude that Ext T(A) is a connected subset of V \ d\(A). But this implies that 
ExtT(A) C V(A) or ExtT(A) C V \ V(A), as claimed. 

(iii) This follows immediately from (i) and the third statement of Lemma 15.11 □ 
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5.4 Isoperimetric Estimates 

We need the notion of diameter, dianry, of a contour 7. To this end, we consider sets Sj^ of 
the form 

= {xe V: x l = k}, 

and define = ^(7) as 

Ii = {k G Z/LZ: n 7 ^ 0}. 

By the fact that 7 is a connected subset of V, we have that Ii is a set of consecutive integers 
mod L (i.e., it is a connected subset of Z/LZ). We define the diameter, diam,7 of 7 in the 
direction % as the number of points in Ii, and the diameter of 7 as diam7 = max i=lj ^ diam^. 

Lemma 5.7 For every contour^, we have 

||7|| > 2 diam7 (5-5) 

and 

|Int7 fl V\ < -H7II diamy . (5.6) 
2 

Proof. Let = Sj^ fl V, and let be the set of edges xy G E such that both x and y 
lie in Sj?\ Consider the configuration A which has 7 as its only contour. If A fl E^ = E^\ 
then Sj. C V(A) (this follows immediately from the fattening procedure used to define 
V{A), and if A n £^ = 0, then C V \ V(A) (this follows from Lemma EH]). In either 
case, 7 must have an empty intersection with Sjjf . (In fact, 7 must have distance at least 
1/4 from this set). For k G Ii, the set £7^ therefore must contain at least one edge in 5 A, 
implying that 7 has at least one intersection with the edges in E® . In fact, by a simple 
parity argument, there must be at least two such intersections. This immediately implies 
that || 7 1| > 2 diamj(7) for all i = 1, . . . , d, which proves the bound (15. 5p . 
Consider now the sets Vi — {x G V | Xi G Ii} and 

Vi = {x G V: 3k G h s.t. \xi -k\< 3/4}. 

Since Ii consists of consecutive integers mod L, the set Vi is a connected subset of the discrete 
torus V, and Vj is a connected subset of V. If dianij7 > 0, the set Vj is non-empty, and 
7 C Vj. If dianij7 = 0, all the edges intersecting 7 must be in the direction i, and since 7 is 
connected, they must all lie in one plane of the torus. In other words, there must be a set 
of the form {x G V: \x{ — k + 1/2| < 1/4} such that 7 is contained in this set. We denote it 
again by Vj . 

While it is in general not true that Int 7 C Vj , we claim that this is true if there exists 
an interface S that is compatible with 7. Indeed, let S be such an interface, and let Si be 
one of the sets S^. C V \ Vj (if there is no such set, Vj = V and there is nothing to prove). 
Repeating the proof of Lemma I4.3[ we see that Si and S must lie in the same component of 
V\7, which by Definition H3J must be the exterior of 7. This proves that Ext 7fl(V\Vj) 7^ 0. 
Taking into account the connectedness of 7, this proves that Int 7 C Vj. 
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For contours whose exterior is defined by the existence of a compatible interface S, this 
immediately implies the bound (15. 6p . Indeed, assume that there exists an interface S such 
that S C Ext 7. Without loss of generality, let us assume that the first component of the 
winding vector N(5) is 1. This implies that every line in the 1-direction intersects S at least 
once. Hence, for all x G W = Int 7 n V, a line through x in the 1-direction intersects 7 at 
least twice. Since W C V 1; this shows that 

|Int7nV A | < -||7|| diam.17 , 

which implies (15. 6p . 

We are thus left with contours 7 for which \ W\ = | Int 711 V\ < L 2 /2. For these contours, 
the isoperimetric inequality of Bollobas-Leader [B] implies that 

INI > \d cdgc W\ > min 2i\W\ 1 - 1/i L d/i - 1 > 2min|W| 1 ~ 1/i L' i/i ~ 1 = 2\W\ 1 - 1/d , 

i=l,...,d i 

where we used the notation d ed g C W for the edge-boundary of W. To complete the proof of 
(15.61) . we will want to show that 

\W\ = |Int 7 n V\ < (diam7) d . (5.7) 

This bound is trivial if dianry = L, so let us assume that dianry < L. Let C = f] i=1 d^i- 
Since 7 C C and V \ C is connected, we must have that either Int 7 C C or Ext 7 C C. In 
the first case, the bound (15.71) is again trivial, and in the second case we use that 

|Int 7 nU| < |Ext 7 nl/| < |Cn V\ < (diani7) d , 

in completing the proof. 

In addition to the above lemma, we will also need bounds on the number of contours and 
interfaces of a given size. More precisely, we need the following lemma. 

Lemma 5.8 There exists an absolute constant C < 00, such that the following statements 
hold 

i) Let 70 be a contour or an interface, and let k > 2. Then the number of contours 7 
such that ||7|| = k and dist(7, 70) < 1/2 is at most \\ r y \\(Cd) k . 

ii) Fix k > L d ~ l . Then the number of interfaces S in V such that \\S\\ = k is at most 
L(Cd) k . 

Hi) Fix a vertex x G V. Then the number of contours 7 such that ||7|| = k and x G Int 7 
is at most (Cd) k . 

Proof, i) Clearly, there is a 1-1 correspondence between contours 7 and the set of edges 
E 1 intersecting them (taking an edge twice if 7 has two intersections with it, so that ||7|| 
is equal to the number of edges in E^), with an analogous statement holding for interfaces. 
Defining a suitable neighborhood relation on the edges in El^, the set E 1 is a connected 
subset in a graph of maximal degree bounded by Kd for some K < 00, and E 1 U E l0 is 
connected if dist(7, 70) < 1/2. Using standard results on the number of connected sets in a 
graph of given maximal degree, we immediately obtain statement i). 
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ii) Let Xq be an arbitrary point in Vl,cLi ^ t> e the straight line through the Xq in i 
direction, and let Ei be the set of edges whose line-segments lie in Lj. If S has non-zero 
winding number in the direction i, then S must intersect each line in direction % at least 
once, implying that Eg must contain at least one edge in E^. Taking Sq to be the union 
Ei U- ■ -UEd, we see that an arbitrary interface S must contain at least one edge in Sq. Since 
So contains dL edges, the result ii) now follows by the argument used in the proof of i). 

iii) Let x = (x±, . . . , x^), and for r = 1, 2, L, let W r be the cube of all points y G V 
such that — | < i/i — Xi < | for all i — 1, . . . , d. Choose R to be the largest r such that 
W r C Int7. Then 7 must intersect one of the 2dR d ~ 1 edges joining Wr to Wr + i, implying 
that the number of contours in question (corresponding to a fixed R) is at most 

2dR d -\Cd) k ~ l < 2d\W R \(Cd) k - 1 < ^k 2 (Cd) k -\ 

where in the last step we used the previous lemma to bound \Wr\ < ^||7|| 2 = To complete 
the proof, we will have to sum over R. But R d = \Wr\ < ^ implies that R < k; thus the 
summation over R gives another factor of k, leading to the bound 

^(Cd)^ 1 < (8Cd) k 
for the number of contours 7 such that ||7|| = k and x G Int7. This proves iii). 

5.5 Matching Contours and Interfaces 

In this section we will show that the partition function Z can be written as a sum over sets 
of matching contours and interfaces. To this end, we establish a sequence of lemmas. We 
start with the following lemma. 

Lemma 5.9 Let 7 be a contour. Then there exists a configuration A such that S(A) = 
andT(A) = {7}. 

Proof. By definition, there exists a configuration A 1 such that 7 is one of the contours 
corresponding to A\. Thus 7 is a connected component of <9V (Ai), and hence a connected 
component of the boundary of one of the components, C, of \(A\). Let A2 be the set of 
edges in A\ whose endpoints both lie in C. Then (V(A 2 ),A 2 ) is a component of (V,Ai), 
and, by Lemma [5.11 (iii), the set V(A2) is nothing but the component C of Y(Ai). 

Thus 7 is a component of d\f(A2). Consider now the complement D = V \ V(A2), and 
its components D , . . . D fc . Then 7 is the boundary of one of these components, say D . Let 
A Q be the set of edges whose midpoint lies in V \ D . By Lemma [5.11 (iv) and (v), we have 
that D = V \ V(A ), which in turn implies that 7 = d\(A ). □ 

Lemma 5.10 Let S be an interface network. Then there exists a configuration A such that 
S(A) = S and T(A) = 0. 

Proof. The proof is identical to that of the previous lemma. 

Recall that each contour has an ordered and a disordered side. We call 7 a contour with 
external label ord if the side facing Ext 7 is ordered. Otherwise it is called a contour with 
external label dis, for disordered. 



24 



Lemma 5.11 Let A G Q. If 7 is a contour with dist(V(A),Int7) > 1/2 and external label 
dis, or a contour with dist(V \ V(A),Int7) > 1/2 and external label ord, then there exists 
a configuration A' C Q such that T(A') = T(A) U {7} and S(A) = S(A'). 

Proof. Consider first the case that i = dis. Define A\ to be the set of edges with both 
endpoints in Int7. By Lemma 15.11 (iii), we have that \{A U A\) = V(A) U \(Ai) = 
V(A) U Int7, which shows that iU^i is the desired configuration. 

For i = ord, we define E\ to be the set of edges whose midpoint lies in Int7. Using 
Lemma [5. II (iv) and (v), we now conclude that A\E\ is the desired configuration. □ 

Corollary 5.12 Let A C Q. Then the contours and interfaces corresponding to a con- 
figuration A G Q are matching. Conversely, any set of matching contours and interfaces 
corresponds to exactly one configuration A £ Q. 

Proof. The first statement is obvious. The second follows from Lemmas 15.101 and the 
previous lemma by induction on the number of contours. Indeed, by Lemma 15.51 the partial 
order on contours leads to a forest on any set T = {71, ...,7 n } of pairwise compatible 
contours, in such a way that 7' < 7 whenever 7' is a child of 7. Adding the interface network 
as a common root, we may proceed by induction from this root to add contours in such a 
way that the new contour added always obeys the condition from the previous lemma. □ 

Note that this corollary, together with the representation ( 14. 5 p for the weights of a con- 
figuration A, established the representation (I4.6P for the partition function Z. 

6 Key Ingredients for the Lower Bound 

As explained in Section 12.31 we will prove our lower bound on the mixing time by proving an 
upper bound on the conductance. This in turn will require both a bound on the probability 
of the set of configurations with at least one interface, and a large deviation bound on 
configurations for which the joint exterior of all contours contains less than (1 — a)L d points 
in V, see Lemmas 16.11 and 16.21 below. We will start with the decomposition of the partition 
function into terms with and without interfaces. 

6.1 Decomposition of the Partition Function 

Let 

funnel = {A E £1 \ S(A) ^ 0}, (6.1) 

On-d = {A E n \ n tunncl | Ext T(A) C V(A)} (6.2) 

and 

n dis = {A E n \ funnel I Ext T(A) C V \ V(A)}. (6.3) 

By Lemma 14.51 (ii), Q = fl ov d U fidis U O tunne i. As a consequence, the partition function 
Z = J2Aen w (A) ( see (14.61) ) can be decomposed as 

Z = Zdis + qZ or <i + ^tunnel) (6-4) 
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where 



J ord 

q 



- ^ W ( A )i Z dis = W ( A ) alld funnel = ( 6 " 5 ) 



Aen OTd A&n dis Aen tu 



nnel 



Note the extra factor of q in (16.41) . which accounts for the fact that there are q different 
ordered phases. 

The results of this section are summarized in the next two lemmas. The first is a finite- 
size scaling bound analogous to those proved in 0, [311 EI] • The second is a large deviations 
bound. 

Lemma 6.1 For all d > 2, there are constants c > 0, qo < oo and Lq < oo such that the 
following statements hold for q > qo and L > Lq : 
(a) If/3 > fa, then 

K^tunnel) < (6-6) 

and 

KO,rd) > — 7— r — e" c/3i . (6.7) 
q+ 1 



(b) J//3 = p , then 

v(Q^a) — 

q + l 



v(n ovd )-^— <e" c/3L (6.J 



To state the next lemma, we define 

_ r /i r- n . . Ixp^-+ r/^nvl s n \ r d 



= {Ae n ord : |Extr(A) n v\ > (1 - a)L d }, 

fi^ = {A G fi dis : |Extr(A) nV\>(l- a)L d }. 



(6.9) 



Lemma 6.2 Let d > 2 and < a < 1. T/ien t/iere are constants c = c(a) > and 
g = 9o( tt ) s,ucft that for q > q and ft > /? we have 

u(Q 0ld \^l) <e- c ^~ X (6.10) 

and 

K^di S \^2) <e-^ Ld -\ (6.11) 



In order to prove Lemma [6.11 we will need upper bounds on Z tunnc i and Z d [ S , as well as 
upper and lower bounds on Z ord . Since contours and interfaces are suppressed if /3 (and hence 
k) is large, the leading configurations to Z ord and Z dis are those without contours, giving a 
contribution of e~ e ° rdL and e~ £disL , respectively. For Z tumie i, the leading configurations have 
a single pair of parallel interfaces of area L d ~ x each, and no contours, giving a contribution 
of at most c~ 2kL max{ge -e ° rdL , e~ £disL }. If we took only the leading configurations into 
account, we therefore would get that Z tunnc \/Z is exponentially suppressed like e~ 2KL , as 
required for the first bound in Lemma 16.11 But of course, this is too naive, since subleading 
contributions have to be taken into account. A systematic way to do this is provided by the 
powerful theory of Pirogov and Sinai. 
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6.2 Ingredients from Pirogov-Sinai Theory 



In this section we will prove Lemma 16.11 To this end, we will express Z tunnc \ in terms of 
partition functions that are analogs of Z ord and Z d \ s for a subset A C V such that 

A is a connected component of V \ dV(A ) for some A £ Q. (6-12) 

We say that a contour 7 is a contour in A if dist(V(7), V \ A) > 1/2, and we say that a set 
of contours T with matching labels has external label £ e {ord, dis}, if the external contours 
in T have external label £. We then set 

Z dis (A) = g c ( v - d ) e - edis|VdisnynA| e- e - d|v - dnynA| Y[ e~ KM , (6.13) 
r 7 er 

where the sum goes over sets of contours V with matching labels such that the external label 
of T is dis, and similarly for Z ord (A): 

qZ OTd (A) = q< v ^e- e ^ v ^ nVnAl e- e ^ lv ^ nVnAl JJ e~ K ^. (6.14) 
r 7 er 

Note that these partition functions are indeed generalizations of the partition functions Z d { s 
and Z ord introduced in (16. 5ft . To see this, it is enough to compare the weights in (16. 131) and 
(I6.14p to the weight w(A) from (14. 5p . which shows that Z dis (V) = Z dis and Z ord (V) = Z ord . 
With the above definitions, the partition function Z tunnc \ can be rewritten as 

funnel = ^ J] e " K " S " II H (^ord(A)) , (6.15) 

S S£S AeCdia(S) AGC ord (5) 

where the first sum goes over interface networks, while C ovd (S) is the set of components 
of V \ Uses S w ith £(C) = ord, and C dis (S) is the set of components of V \ Uses ^ with 
£(C) = dis. The formal proof of this, by now almost obvious, identity uses again the forest 
structure of sets of pairwise compatible contours established in Lemma [531 and is similar to 
that of Corollary [Ell 

6.2.1 An alternative representation 

We will need a representation for the partition functions Z dis (A) and Z ord (A) which does 
not involve the restriction of matching labels. To this end, first sum all terms in (16.131) and 
(I6.14p which lead to the same set, r ext , of external contours. Taking into account the forest 
structure of sets of pairwise compatible contours established in Lemma 15.51 this leads to the 
identities 

Z ord (A) = ^e^ e - d l AnExtr ^ ny l J] e- K "^Z dis (Int 7 ) (6.16) 

Text 7£F e xt 

and 

Z dis (A) = ^e- e -l AnExtr - tny l J] e-*MgZ ord (Int 7 ), (6.17) 

Text 7£F ex t 
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where the sums run over sets of mutually external contours in A which all have external 
label ord and dis, respectively. 
Defining 

K ord 7 = e " ll7 " „ : T 7 \ , 6.18 
Z ord (Int7j 

we rewrite ( 16 . 1 6[) as 

Z ord (A) = ^e- e ^l AnExtr - nV/ l J] K ord (7)Z ord (Int 7 ). (6.19) 

Text 7€r ex t 

Inserting (16. 191) inductively into itself, and using the forest structure of sets of pairwise 
compatible contours one last time, we finally arrive at the representation 

-Z'ord(A) = e ~^\Anv\ J- Y[ K ovA (l) , (6-20) 

r 7 er 

where the sum now runs over sets T of pairwise compatible contours in A which all have 
external label ord. In a similar way, one shows that 

Z dis (A) = e - e ^ Any \ II ( 6 - 21 ) 

r 7 er 

where the sum runs over sets V of pairwise compatible contours in A which all have external 
label dis, and 

ifdb(7)=e -*«^a^). (6 . 22) 

Z dis (Int7) 

The representations ( 16.201) and ( 16.211) give Z ord (A) and Z d ; s (A) as partition functions of 
the so-called abstract polymer systems, see, e.g., [13] or [7], for a review. As a consequence, 
the logarithms of Z ord (A) and Z d i S (A) can be analyzed by absolutely convergent expansions 
(so-called Mayer-expansions), provided the weights (16.181) and (16.221) are sufficiently small. 

The following lemma gives the bounds needed to apply these expansions. Given the 
geometric preparations of the last section, its proof follows from a careful extension of the 
methods of PU E]. For the convenience of the reader, we give it in the appendix. 

Lemma 6.3 Let d > 2. Then there are constants qo > and c > 0, as well as two real- 
valued functions / ord = /ord(<?>/3) and fd is = /dis(^, /?) such that the following statements hold 
for £ G {ord, dis} ; q > q , and (3 > (3 : 

(i) Let f = min{/ ord , / d i S }, and let ag = fg — f . If 7 is a contour with external label £ and 
a^dianry < cf3, then 

Kiin) < e~ M . (6.23) 

(ii) If A C V is of the form (16. 12[) . then 

Z e (A) > e -(h+^)\AnV\ e -\\dA\\ (6>24) 

and 

Zt(A) < e (-/+^)|AnV| e 2||9A|| maxe -^|AnExtr cxt ny| " e - f /j M> (g 25) 

Text _ 

7er cx t 
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where the maximum goes over sets of mutually external contours in A which all have external 
label £, and e L = 2e~ c/3L . 

(Hi) ford < /dis if P > 0o, with equality if f3 = (3 . 

6.2.2 Proof of Lemma EH] 

We start by noting that by (11. 5ft and ( 14.41) . the assumption /3 > /3q implies that 

K > > — : log q , 

~ 2 4 - 2d &H 2 ' 

provided q is large enough (depending on d). We also recall the notation V = = (7LjUE) d 
for the vertex set of the d-dimensional discrete torus of sidelength L. 
To prove (16. BJ) , we combine (I6.15P with (I6.25p . to conclude that 



funnel <J2H e (~ f+e L )\^V\ e 2\\dX\\ ge (-/ +££ )|AnV| e 2|| S A|| 

S SeS AGC dis (5) AGC ord (5) 

s ses 



(6.26) 



where the sum goes over interface networks. In the last step we used that each interface 
bounds precisely one ordered and one disordered component. Using the facts that q < e 2dK e d , 
that there are at most 2dL{Cd) k interfaces of size k in Vj,^ (for an appropriate universal 
constant C < oo), and that the sum over interface networks contains at least two interfaces, 
the bound (16.261) implies that for q and L large enough (depending on d), we have 



Wnei <e^- f+ ^ Ld J2( E e 2dK+d 2dL(Cd) k e- { 

n>2 fc>Ld-i 



rd-X v n 

e (-f+e L) L« y | e 2 M 2dL ^— A , (6.27) 




e -l KL d~l n < e _ fL d e _l L d-l ^ 



n>2 

where we used that L d ei < e^ 2d ~ c ^ L to control the factor e £Lljd . 
Applying the bound (I6.25P to Z dis = Z^ s {y Ljd ), we get 

^dis < e-f Ld+Ld ^ < e^ L \l + e~ c ^ 2 ) , 

while the bound (I6.24p . together with the fact that / or d = / if /3 > j3o, gives 

Z OTd > e~f Ld ~ Ld ^ > e^ Ld (l - e- c ^ 2 ) . 



Together with ( 16.41) . this gives the first statement of the lemma. 

To prove the second statement, we use that / = /dis = /ord if /3 = /3 which, together 
with (E2U) and flQfl gives 



logZ dis + fL d \<L d e L <e- clSL/2 and | log Z ord + fL d | < L d e L < e 



■c/3L/2 
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6.3 Large Deviation Bounds 

We now prove Lemma [6.21 by starting with the proof of (16.101) . Let 



« = \ E w ^ 



Q 

\ExtF(A)nV\<(l-a)L d 

so that 



Proceeding as in the derivation of (16.201) . we rewrite 



K^ord \ = ^ < Si . (6 .28) 



7<ct _ -e ord L' 
J ord c 



e n j ^(7) 

T: |ExtrW|<(l-a)L d 76r 



where the sum runs over sets T of pairwise compatible contours in Vl^, all of which have 
external label ord. 

Let h be an arbitrary non-negative number. Using the fact that |Ext r D V\ < (1 — a)L d 
implies that X] 7 er |I n "t T H ^| > aL d , we then bound 

Z<« < e~^ Ld £ ^ |Int ^- Qid) J] K Ui) = e- ahLd Z^ d , (6.29) 
r 7 er 

where 

^ = e~ e - Ld Eri e/l|Int7ny| ^(7). 

r 7 er 

Next we estimate the dependence of Z^\ on h. To this end, we set lift (7) = e h ' Int7ny ' K Q rd{l) 
and use a Peierls type argument to bound the derivative of Z^\. Explicitly, we first rewrite 
the derivative of Z^\ as 



^ QT |int 7 nv|)n ^(7) 



— z 

r 7 er 7 er 



-e ord L d 



E |Int 7 nV|^n K h {i) 
jcv L:d rs 77 'er 

= e -e ^ j2 |int T ny|^( 7 )E / II K »M, 

7 cVi, d r' 7 'er' 

where the sum over 7 C Vl,<i denotes a sum over contours in Vi,d with external label ord, 
and the sum ^ ' denotes a sum over sets V of pairwise compatible contours in Vl,cl such that 

• all contours 7' G V have external label ord; 

• all contours 7' G V are compatible with 7. 
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Removing the second constraint, we bound this sum by 

e' n k »m * e n k *w = ^ dLd ^i 

v 7'er' r 7 'er 

which shows that 

d . _cm 1 d 



^ord 7CVt 



| 7 ||2 e fc||7||L e -e)9||7|| 

where we used the bounds (15. 5p . (I5.6p and (I6.23p . together with the fact that / or d = / for 
j3 > f3 , in the last step. 

Now assume that h < c/3/(2L). Since there are most dL d (Cd) k contours of size k in Vl j( 2, 
we conclude that 

l l0 S^d < \52dL*(Cd) k k 2 e-°W 2 < \L\ 

k>2 

and thus 

y{h) 7 ahL d /4 
^ or d S ^ord 6 , 

provided /3 is large enough (depending on d and a.) Inserted into ( 16.29j) and ( 16.28[) . this 
gives 

K^d \ ^ rd ) < e" 3 ^ 4 = e-^-Vs j 

where we have set h equal to cf3/ (2L) in the last step. 

This proves (16.101) . For a dis L < c/3, the bound (16. lip is proved in exactly the same way, 
but for a dis L > eft, this strategy does not work, since (16.231) is not at our disposal anymore. 
In stead, we use (16.241) . (I6.25P and (15. 6p to bound 



■^ord r<!xt ^ 



< e 2L ^ maxe-^l Extr - ny ' ' e ~^ Int ^ 

Text „ 

7er cxt 

< e 2Ld ^e-£ Ld , 



where we used that |Ext r ext Pi V \ = L d — X] 7 er cxt \^ n ^ 7 ^ ^1 m the last step. This concludes 
the proof of the lemma. 



7 Lower bounds on the mixing time 

Lemma [6.11 and Lemma [6.21 give us the necessary ingredients to prove our main result. We 
start by proving a key theorem which expresses the statements of these lemmas in terms of 
the probability measure /x on spin configurations, rather than the random cluster measure 
v. 
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7.1 An Important Probabilistic Estimate 

Recall that for a spin configuration cr e [q] , we defined -E(cr) to be the set of all edges 
whose two endpoints have the same color. We also introduce the set, C(cr), of connected 
components of the graph (V,E(cr)), that is, the set of monochromatic components. 

Theorem 7.1 Let d > 2 and < a < 1/2. T/ien there are constants c > 0, qo < oo and 
L < oo sitc/i for q > qo and L > L the following statements hold: 
(a) If (3 = (3 , then 

JadL d <\E{a)\<{l - a)dL d ^ < e"^^ 1 (7.1) 



and 



(b) If/3 > fa, then 



and 



/i(|£(<x)| > (1 - a)dL d ) - < e~^ L . (7.2) 



H[ aL d < max \V(C)\ < (1 - a)L d ) < e~ cfiL (7.3) 



( max > (1 - «)L d ) > - e~^ L . (7.4) 

\C£C(<t) J q + 1 



Proof. To relate the statements of Lemmas 16.11 and 16.21 to the theorem, we use that both 
the spin measure /x and the FK-measure v are marginals of the Edwards-Sokal measure ir. 
Consider thus a configuration (cr, A) with positive measure vr((cr, A)). Under this condition, 
all spins in a component of (V, A) must have the same color, implying in particular that 

A C E{a) . (7.5) 

It turns out that with high probability, \A\ is not much smaller than |-E(cr)| either. More 
precisely, we will prove that 

7r(|£(cr)| > \A\ + adL^j < e ~( & ^) dLd for all a G (0, 1). (7.6) 
We will also show that, again under the condition that vr((cr, A)) > 0, 



and 



max \V(C)\ > max \V(C)\ (7.7) 

CeC(o-) C£C(V,A) 

max IV (C) I < max \V(C)\ + \E(tr) \ A\ + 1. (7.8) 
cec(cr) ceC(v,A) 

As we will see below, these bounds, together with Lemmas 16.11 and 16.21 imply the statements 
of the theorem. 

Before showing this, we will prove (17.61) - (17. 8p . We start with the proof of (17. 6p . To this 
end, we rewrite the left hand side as 

7r(|£(er)| > \A\+adL d ^j = n(\E(tr)\ > \A\+adL d <r)/x(<r). 
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But given cr, a configuration A according to the conditional measure 7r(- | cr) is obtained 
by deleting the edges in E(er) independently with probability e _/3 . The number of deleted 
edges is therefore equal to a binomial random variable with parameters m and e _/3 , where 
m = \E(cr)\ < dL d . We now bound the probability that the number of deleted edges X is 
larger than adL d as follows: 

Pr(X >adL d )= J2 (™)e-P k (l-e-P) m - k <e-P &dLd 2 dLd . (7.9) 

k>adL d ^ ' 

This implies the bound (17. 6p . 

Next, we observe that (17.71) follows from (17. 5p . To prove (17. 8p . we first show that any for 
C G C(a) and any D C E(C), we have: 

\D\-\V(D)\<\E(C)\-\V(C)\ + 1. (7.10) 

The simplest (albeit non-elementary) way to see this is by recalling that the dimension of 
the cycle subspace of the edge space of a graph on m edges and n vertices and k components 
equals m — n + k; we view here the collection of all edges as the m-dimensional vector space 
over GF(2) and consider the edge sets of all simple cycles as a subspace of it. For non-empty 
D, this gives in fact the stronger bound 

\D\-\V(D)\ + 1<\E(C)\-\V(C)\ + 1. 

If D = 0, the bound (17.101) is trivial, since the right hand side of (I7.10P is non-negative by 
the fact that C is connected. 

Consider now a component C of (V(A), A). Since A C E(cr), there must be a component 
C E C(tr) such that V{C) C V(C) and E(C) C E(C). Applying (17301) with D = E(C), we 

^ \V(C)\<\V(C)\ + \E(C)\-\E(C)\ + 1 

< \V(6)\ + \E{<r)\-\A\ + l 

< max \V(C')\ + \E(tr)\ - \A\ + 1. 

CeC(v,A) 

This in turn implies (17. 8p . 

We are now ready to prove statement (a). To this end, consider a configuration A £ 

where (and fi^Ts) are as defined in ( 16.91) . Then ExtT(A) D is connected set by 
Lemma [4.51 If cr is such that ir((cr,A)) > 0, all edges joining two points in ExtT(A) D V 
must then be part of E(cr). Since the number of edges intersecting the complement of 
Ext T(A)DV is at most 2d(L d - |Ext T(A)nV\) < 2daL d , we concluded that E(o) contains 
at least dL d — 2daL d edges. In summary 

AeQ^ and vr((cr,A)) > |£?(<r)| > (1 -2a)dL d . (7.11) 

On the other hand, by the fact that |A| < d|V(A)| for all A C E, we have that 

AeQffi =^ \V\ V(A)\ > (1 -a)L d =^ \A\<adL d . (7.12) 
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We now turn to the first bound of the theorem. To this end, we use (17. lip . (17. 6p . and 
( ELTgp to bound the left hand side of dZHJ by 

K^tunnei) + 77 (a G ^ord and |£(<x)| < (l - a)dL d ^j + 7r(A G fidis and |£(er)| > a<iL d ) 

< ^tunnel) + K^ord \ ^f) + K^iis \ >) + e^^ 1 )^. 

Bounding the terms on the right hand side with the help of (16. 6p . ( I6.10P and (16. lip , we 
therefore obtain that there exists a constant c > depending on a and d such that 

fi(adL d < \E(tr)\ < (1 - < 4e" 2c/3L " 1 < e"^" 1 , 

provided q (and hence /3) is large enough. This proves the bound (17.11) . 

The proof of the bound (17.21) is similar. Indeed, starting again with the implication 
(17. lip , we have 

n(\E(<r)\ > (1 - a)dL d ) > v(n^) > _ e-^" 1 - 

where we used the bounds (16.71) and (16. 10[) in the last step. 
On the other hand, by (17.61) . 

n(\E(tr)\ > (1 - a)dL d ^j < u(\A\ > (1 - 2a)dL A j + e ~ {a ^ 1)dLd 

< funnel) + v(i&) + K^dis \ fig. - * ) + e'^^ . 

Combined with the bounds (16.61) . (16. 8p . (16. 10[) and (16.111) . this provides a matching upper 
bound on fi(\E(cr)\ > (1 — a)dL d ^j, completing the proof of (17. 2 j) . and hence of part (a). 

The bounds of part (b) are proved in a similar way. Indeed, let A G and let 

A ext be the set of edges with both endpoints in ExtT(A). By Lemma [4.51 (hi), the graph 
(V(A ext ), A ext ) is connected, implying that maxceC(v,A) \V(C)\ > (1 — a)L d whenever A G 
^ird- Taking into account the bound (17. 7[) we get that 

A G Q^l and it ({a, A)) > => max \V{C)\ > (1 - a)L d . (7.13) 

CEC(cr) 

Together with (16 . 7|) and (I6.10p . this immediately gives the bound (I7.4p . 

We are thus left with the proof of the bound (17. 3p . To this end, we note that if |^(A)| < 
\aL d and |i?(cr) | < | + \aL d —\ then, by (1 7 . 8 1) . the largest component of (V,E(cr)) has 
number of vertices less than | V^(^4) | + |-E(cr)| — \A\ + 1, which in turn is less than aL d . 
Combined with the bound (17. 6 j) . we conclude that 

7r (a G {l { */ 2) and max \V(C)\ > aL d ) < e -(^-^ d 
V ms c*ec(<r) / 

where a = ^ — -r^. Combined with ( I7.13p . this implies that the left hand side of (17. 3p can 
be bounded by 

,i(fl\(n2u^)) + e-^. 

Together with the bounds (16. 6p . (16.101) and (16.111) . this gives the desired bound (I7.3p . □ 
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7.2 Proof of the SW bound in Theorem 11.21 

Let (3 = (3q. Let S = {cr : |-E(cr)| > (1 — a)dL d }. We will show that $s is exponentially small 
in /3L d_1 , which will establish the theorem. For q (and hence j3o) large enough, tt(S) > 1/2, 
using (17. 2p . Also, n(S c ) > 1/q — e _c/3L > l/2q, if L is large enough. Thus 

*^M$k i4qQis - s ' r (7 - i4) 

Let S = {cr : adL d <\E(cr)\<{l - a)dL d }. Then 

Q(S, S c ) = Q(S, S c \ S ) + Q(S, So). (7.15) 

Now 

Q(S,S ) = Q(S ,S) < ti(So) < e- cf3Ld -\ 

using (17. ip . while 

Q{S,S C \S ) =tt{S)Pi(\E((t')\ <adL d \\E((r)\ > (l-a)dL d ), 

where cr is a /i-random spin configuration, and cr' is constructed from one step of the SW 
algorithm. The above probability is in turn the probability that at least (1 — 2a)dL d edges 

are deleted in one step of the SW algorithm, which is at most {2e^^^ 2a ^ , using (17. 9p . 
Choosing a = 1/3, and q large (so that (3 is large) and L sufficiently large, yields: 

Q(S, S c \ S ) < tt(S) [2e~^- 2a) ) dLd < e~ c ^ Ld ~\ 
This implies the desired bound on conductance, 

$(P SW ) < e- c ? Ld -\ 

which together with (12.91) concludes the proof of the lower bound (ll.lOp on the mixing time 
of the SW algorithm. □ 

7.3 Proof of the HB bound in Theorem 11.21 

Let cr G [q] v and let C be a component of (V, E(cr)). We say that C has color k, if a x — k 
for all k G V(C), and denote the set of components of color k by Cfc(cr). For k = 1, . . . , q 
and < a < 1/2, we then define 

h ( k a) = |cr G [q} V : 3 C G C k {a) s.t. \V{C)\ > (1 - a)\V\). 
Note that the sets Vt^ are mutually disjoint, so by symmetry and the bound (17.41) . we have 

^)>-L_-le~^. (7.16) 
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Finally, let 

= {* e [q\ v : T^x ) \V(C)\<o l L d }. 

We complete our proof by estimating $5 (see (12.71) ) for S = First notice that for q 

(and hence also (3) sufficiently large fi(S)fi(S c ) > l/4q so that 

® S <MQ(S,S C ). 

Since the heat bath algorithm can only change one vertex at a time, it does not make 
transitions between the different sets Q k . For a small enough, it cannot make transitions 
between and f^ s either. Indeed, changing the color of a single vertex can not break 
a component C G C(cr) into more than 2d + 1 new components: in the worst case, C gets 
broken into a single component of size one and 2c? components of size ((1 — a)L d — 1)/ (2d). 
For a sufficiently small (say a = l/(4d)), the heat bath algorithm therefore cannot make 
transitions between and . Defining So as the set of configurations which are neither 
in Vl^ nor in one of the sets fij^, we thus have 

Q(S, S c ) = Q(S, S ) = Q(S , S) < n(S Q ) < 

where we have used the bound (17.31) in the last step. Recalling that [3 > (3 = ^logg + 
0(q~ l / d ), we see that for L sufficiently large, we have $5 < AqQ(S, S c ) < e~ cl3Ld 1 ^ 2 , as 
required. □ 
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A Technical Estimates using Truncation 

Throughout this section, we will assume that 

> maxjd log(dC), ^logg - l}, (A.l) 

where C is the constant from Lemma 15.81 and C\ is a suitable constant to be chosen in the 
course of the proof. In fact, we will prove statements (i) and (ii) of Lemma I6TB1 for all such 
that flA~Tj) holds (whether /3 > /3 or not). 

Also, for the purpose of this appendix, we will use the symbol |A| for the cardinality of 
the set Vfl A, so that expressions of the form | VnExt Tfl A| can be simplified to | An Ext T|. 

A.l Truncated Contour Models 

We need some preparation. We start by bounding the factor ge _K " 7 " appearing in the weight 
( I6.22p of a contour with external label dis. To this end, we observe that the smallest contour 
with disordered external label has size ||7|| > Ad — 2. Combined with the assumption (lA.ip 
and the assumption d > 2, this gives 

g e -«|[7ll < e ^+i) e -«[| 7 || < exp ((^±| - «)|| 7 ||) < e -flMI , (A.2) 

where in the last step we assume that C\ is chosen large enough to guarantee that k = 
§ + O(e-0>i + iM. 
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As usual in Pirogov-Sinai theory, we next introduce a truncated model. It is given in 
terms of the truncated activities 



Kfa) = min{^( 7 ),e-^- c/3+1)l1711 } 
and the corresponding partition functions 

Z;(A) = e-^ A l^n^(7), 

r 7 er 



(A.3) 



(A.4) 



where I G {ord, dis}, the sum in (1A.4|) runs over sets T of pairwise compatible contours in A 
which all have external label £, and c is a small enough constant; we will choose c = 1/20, 
implying in particular that K' g (j) < e _c ^" 7 ". 

Let x G V, let 70 be a contour or an interface, let i G {ord, dis}. With the help of 



Lemma [5.81 we then bound 



£ ^( 7 ) e c^)iwi< £ 

7:Int79z 7:Int79z 



e -(|-2 C /3)||^|| < 



E 

fc>l 



e v s 



(#-2c/3)fc 



(Cd) fc < 1, 



(A.5) 



and 



-(f-2c/3)fc 



(c^) fc < hoi 



(A.6) 



£ ^( 7)e (^)IMI<|| 7o ||^ e -„ 

7:dist(7,7o)<l/2 k>l 

provided C\ is sufficiently large. The above bounds imply absolute convergence of the cluster 
expansions for abstract polymer systems, which in turn gives the existence of the limits 



ft = lim jf ) with 



-—\o g Z[{V L4 ) 



(A.7) 



where I G {ord, dis} and Vlj denotes the d- dimensional torus of sidelength L, see, e.g., 
[T3| [7] for a review of cluster expansions for abstract polymer systems. These methods also 
imply that, for I G {ord, dis}, and A C V of the form (16.121) . we have \fy — f\ \ < Sl and 



logZ;(A) + /,|A| 



< ||9A|| +e L \A\ 



(A.8) 



where, as before, El = 2e C/3L . We will assume that C\ in (1A. 1|) has been chosen in such a 
way that Lel < 1. 



A. 2 Proof of Lemma 16.31 (i) and (ii) 

We first note that Zi(A) > Z' e (A), so in view of flA.8j) . we have 

Z e (A) > e -(/i+^)|A| e -||9A||_ ( A 9 ^ 

Next we recall that K[(-f) < e" c/311711 . As a consequence ^(7) < e c/3 " 7 " whenever Kg{^/) = 
^(7). To prove Lemma 16.31 (i) and (ii), it is therefore enough to establish the following 
lemma: 
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Lemma A.l Under the condition flA.ljl . we have that 



(i) Ki{pf) = K'fi^f) whenever 7 is a contour with external label i and a^diarrry < c/3. 

(ii) For all A of the form (16.12j) . 

^(A)<e^^l A le 2 " aA llm a xe-^l AnExtr -l J] e"^W , (A 1Q) 

where the maximum goes over sets of mutually external contours in A which all have external 
label I. 

Proof. We prove the lemma by induction on the levels of A and 7: Here the level of a set 
A is defined to be zero if there are no contours 7 such that 7 is a contour in A. The level 
of a set A is defined to be k if the highest level of a contour 7 in A is k — 1, with the level 
of a contour inductively defined to be 1 plus the level of its interior. Note that with this 
definition, the levels of two contours 7, 7' with 7 < 7' differ by at least two. 

Assume that A has level 0. Recall the definition of eg from ( 14. 2 p and (14. 3p . Then 
Zi(A) = e~ e£ ! A ! < e~^' A ' = e~^' A 'e _af ' A ', where the inequality may be seen as follows: first 
for a fixed L, consider Z' e with A being the entire torus of size L d \ the term with no contours 
gives a contribution e~ eiL,d , implying the desired bound for f\ . Taking the limit L — > 00, 
one gets the inequality. This proves ( jA.lOj) for sets A of level 0, the base case. 

Next assume that the bound ( lA.lOp in (ii) has been proven for all sets A of level k or 



less. If 7 has level k + 1 or less and label I = dis, then 



Yy- i \ ■ _£| H iZ ord (Int7) 

^dis(7) < e 8 11711 ) ; 

Z dis (Int7j 

< e (of+2ei,)|Int7| e -(|-3)||7|| <- e -(|-3-|(a f +2e L )diam7)||7|| 



where we used ( 1A.2|) in the first inequality, the inductive assumption ( 1A.10P and the bound 



(IA.9P in the second, and the bound (I5.6P in the last. Bounding e^diam7 by Le^ < 1 
and using the assumption G^diaHry < c/3, this gives K^i^y) < e - ^ -4- ^" 7 " and hence 
Kdisd) = K'disil) ( a g arn provided C\ is sufficiently large). The bound for contours with 
ordered external label is exactly the same. 

Finally, assume that A has level k + 2, that (i) has been proven for all contours of level 
at most k + 1, and that (ii) has been proven for all sets of level at most k. Define a contour 
7 with external label t to be small if (i^dianr/ < c(3, and large otherwise. Consider the 
representation ( I6.13P for Zdi S (A), and fix, for a moment, the set ri argc of all large external 
contours contributing to the right hand side. Summing over the remaining contours, we 
get a factor of Z or d(Int7) for the interior of each contour 7 £ ri argc , as well as a factor 
Z^j™ a11 ^ (Ext ri argc ) for the exterior of ri argc , where Z d j™ all ^(A') is obtained from Zdi S (A') by 
dropping all configurations with large external contours. Thus 

Z dis (A)= ^4 s mall) (Extr largc ) n qe~ K ^Z OTd (Intj), (All) 

-L large 

where the sum goes over sets of mutually external, large contours with disordered external 
label. 
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Since 7 is small whenever 7 < 7' and 7' is a small contour with the same external label 
as 7, the representation ( 16 .2 1 j) for Zj™ a11 ^ (Ext ri arge ) contains only small contours, implying 
that for all these contours i£dis(7) = ^dis(T)- As a consequence, 



4 s mall) (Extr largc )<^ is (Extr 



which allows us to use the estimate (1A.8j) to estimate the factor (Ext ri arge ) in (lA.llj) . 



large J j 

(small) , 



Using the inductive assumption (ii) to bound the factors Z ord (Int7) and the bound (1A.2|) to 
estimate the factors ge~ K " 7 ", this gives 

Z dis (A) < ^ e (eL ~ /dis)|AnExtriar s cl e l|9Extriargcl1 Y\_ e" ( ^" 2)ll7ll e (£i_/)|Int71 

Tlarge T £ -Tl ar g C 

= e ( £i -/)|A| + ||OA|| e -a dis |AnExtr largc | JJ e -(|-3)|| 7 || _ 

^large T^-^largc 

In order to prove statement (ii), we will have to show that 

^2 e -^dis|AnExtr largo | jj e _(|_ 3 _<£)|| 7 || ^ e || 8A || ^ (A.12) 

^large O'^^largc 

To this end, we define 
K{l)~- 

and 



e (f 4 £)h\\ if 7 is a large contour with external label dis 



otherwise, 



r 7 er 

where the sum runs over sets V of pairwise compatible contours in A' which all have external 
label dis. We also define 

f = ~\ogz<y L4 ). 

We will need the following lemma, whose proof we defer to Appendix IA.3I . 
Lemma A. 2 Let A' be of the form (16.121) . Then 

logZ(A') + /|A'| < ||9A'||. (A.13) 

Furthermore —f= \ f\ < 

Recalling the definition of -^(7), we now use Lemma [A. 21 to bound the left hand side of 

$KM b y 
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^ e /|AnExtr cxt | ■ e -||7ll£( 7 ) 



r C xt 7sr ex t 

< ^ e />nExtr cxt | ^( 7 )e/l Int ^Z(Int 7 ), (applying flAH to Int 7) 

Text 7Gr ex t 

= e/|A 'E II ^(7)^(Int 7 ) 

Text 7Sr GX t 

= e /|A| Z(A) < e l|9A|1 , (once again by (EP3| ) , 

proving (IA. 12[) . Note that in the above, the sums run over sets of mutually external contours, 
all of which have external label dis. 

This concludes the proof of (ii) for sets A of level k + 2 and t = dis. The proof of (ii) for 
I = ord is identical. 



A. 3 Proof of Lemma IA.2 



As before, one can use the forest structure of sets of pairwise compatible contours to rewrite 
Z(A') as a sum over sets of mutually external contours in A': 

^ A ') = E II ^(7)^(Int 7 ). 

r cx t 7£r ex t 

For G\ sufficiently large, the partition function Z(A') can again be analyzed by convergent 
Mayer expansion, leading to the bound (IA.13P (the term proportional to El is absent since 
we defined / without taking the limit L — > 00). 
To bound /, we use that 

00 n 

e-f~ Ld = Z(V L , d ) < E^T E = exp£>(7)), 

n=0 ' (7i,...,7„) t=l 7 

where the first sum goes over (not necessarily compatible) sequences of contours in Vl,<i, with 
external label dis. To bound the sum in the exponent, we use that Ki^y) = 0, unless 7 is 
large, which by (15. 5 j) implies that || 7 || > k$ = 2cf3/a&i S . Furthermore, if 7 is a contour in Vlj 
it must be incompatible with one of the L d contours obtained by considering configurations 
with one disordered edge in the 1-direction. Since all of these have size 2, we may use 
Lemma [5.81 to bound the number of contours of size || 7 || = k by 2L d (Cd) k . Using these two 
observations, we get 

£ tf( 7 ) < 2L d £ (Cd) k e-^- 4 ~^ k 

7:||7||=fc>fco 

< 2L d Y,[{Cd)e-^-^] k 



k>ko 



2L d £ [\ e 

k>ko 



< L d e~ c>Sko 
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again provided that C\ is sufficiently large. Thus we have that 



_ f _ I f I < p -c/3k < 1 _ fl dis < ^dis 

1 1/1 " " c/3A; 2(c/3)2 " 2 ' 

A. 4 Proof of Lemma 16.31 (iii) 

We start with the observation that the weights Ke(~f), and hence the weights ^(7) are 
continuous functions of f3. Since the free energies fe are given in terms of an absolutely 
convergent power series in the weights ^(7), they are continuous functions of (3 as well. 
Taking into account this continuity, the following lemma immediately implies Lemma 16.31 
(iii). Recall the definition of M(/3) from the introduction - see below (jl.4p . 



Lemma A. 3 Assume that (lA.lj) holds. 



ft) Ifa ord = 0, then M(/3) > 0. 
(it) Ifa ord > 0, then M((3) = 0. 

Proof. At this point, the proof of Lemma [A.3I is pretty standard. We therefore only sketch 
the main steps. 

First, we note that for A = A^ = {1, . . . , L} d , the representations (12. ip and (12. 2p can 
be generalized to the model with 1-boundary conditions defined in (II. 4p . Indeed, let G + 
be the induced graph on A + = {0, 1, . . . , L + 1} C Z d . The Edwards-Sokal measure 7Ta,i 
corresponding to /iA,i can then be obtained from the measure t\q + by conditioning on a x = 1 
for all x £ A + \ A and xy £ A whenever {x, y} C A + \ A. Next, we observe that in the 
conditional measure ttg + {- | A), a spin at a vertex x £ A has probability 1/q of taking the 
value 1 unless x lies in the same component of (A + , A) as A + \ A, in other words, unless 
x £ Ext (A). Keeping these two observations in mind, the derivation of the representation 
(I6.20p can easily be adapted to obtain a contour representation for the magnetization. Setting 
A = (-1, L + \) d C M d , this gives 



7er 



where the sum goes over sets of pairwise compatible contours in A with external label ord. 
Note that we have chosen A in such a way that all edges in A + \ A lie in Ext 7 whenever 7 
is a contour in A, corresponding to the above conditioning in hg+- 

If a or( j = 0, the weights i^ rd(7) are bounded by e _c/3 " 7 " for all 7. As a consequence, we 
may use a standard Peierls argument to show that the probability that a given point x £ A 
lies not in Ext V is small uniformly in L and x £ A^, implying that M(/3) > 0, which proves 

(i)- 

Assume finally that a or d > (which implies in particular that a dls = 0). We will show 
that with probability tending to one, |A (7 Ext T| < L d ~ e . Since the ratio |A D Ext T|/|A| in 
the definition of the magnetization is bounded uniformly in L, this will show that M\((3) — > 
as L 00. 

Recall the definition of large contours from the last proof, and let ri argc be the set of 
external contours in T which are large. Then |A H Ext T| < |A R Extri arge |, implying that 
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it will be enough to show that the probability that |A n Extri argc | > L d ~ x l 2 goes to zero as 
L — > oo. To bound this probability, we will prove an upper bound on the sum over contours 
with | A H Extri argc | > L d ~ l l 2 and a lower bound on Z or d(A). 

To obtain the desired upper bound, we proceed as in the proof of Lemma IA~T1 (ii). leading 
to the estimate 



e -e ord |A + | J- l[K old (l)< 

< e -/|A+| e 2||aA|| max e -^|AnExtr cxt | TT e -f/3|| 7 || (A.14) 

~~ r cxt : 11 

lAnExtrext^L 1 *- 1 / 2 7€r cxt 

< e -f\A+\ e 2\\dA\\ e -^flL*-i/l 

where the maximum in the second to last line goes over sets of mutually external contours 
in A which all have external label ord. 

To bound Z or d(A) from below we restrict the sum in (16.161) to a single term, the term 
Text = {7o}; where 70 is the contour 70 = <9[l/4, L + 3/4] d . This gives 



Zord(A) > e -eo rd |A + \A| e - K || 70 11^^^) > e -fL* e -2d( K+ l+0(c 



d-l 



■ 



where we used the bound (IA.9D and the fact that adi S = in the second step. After extracting 
the leading contribution e~^ L , the right hand side falls at most like an exponential in L d ~ l , 
while the corresponding decay in (1A.14|) is exponential in L d ~ l l 2 . This proves that the ratio 



of ( 1A.14I) and Z or d(A) goes to zero as L — > oo, as desired, completing the proof of the 
lemma. 
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